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Picosecond  Laser  Interactions  with  Semiconductors 


The  personnel  associated  with  this  project  during  this  reporting  period 
consisted  of: 

Dr.  N.  Bloembergen,  Principal  Investigator 


Dr.  H.  Kurz,  Senior  Research  Fellow 
Dr.  A.  M.  Malvezzi,  Research  Fellow 


for  public  release  and  sale;  its 
distribution  is  unlimited. 


Dr.  D.  von  der  Linde,  Visiting  Scholar  (3/1/84-9/30/84) 


Dr.  G.-Z.  Yang,  Visiting  Scholar  (from  1/1/84) 


Optical  probing  with  picosecond  pulses  at  various  wavelengths  of  the  com¬ 
plex  index  of  refraction  of  ft  in  semiconducting  samples  of  silicon  and  ger¬ 
manium,  in  which  energy  has  been  deposited  by  preceding  picosecond  pulses,  has 
provided  firm  quantitative  experimental  evidence  that  the  energy  absorbed  by 
carriers  gets  transferred  to  the  lattice  on  time  scales  shorter  than  one 
picosecond.  If  the  energy  absorbed  per  unit  surface  area  is  sufficient,  melting 
and  superheating  of  the  fluid  also  occurs  on  such  time  scales.  These  findings 
have  been  published  or  presented  at  international  meetings. 


During  the  past  year  we  have  obtained  evidence  that  a  GaAs  surface  layer 
irradiated  by  a  green  pulse  (X  «  530  nm)  of  20  ps  duration  also  undergoes  a 
transition  to  the  molten  phase  on  a  picosecond  time  scale.  In  this  case  a  dif¬ 
ferent  optical  probing  technique  was  used.  The  incident  light  pulse,  polarized 
along  the  (1,1,1)  crystallographic  direction,  creates  a  second  harmonic 


reflected  ray.  The  second  harmonic  energy  generated  at  X  =  26S  nm  is  propor¬ 
tional  to  the  square  of  the  incident  intensity,  as  long  as  the  sample  remains 
crystalline.  When  a  layer  equal  to  the  absorption  depth  at  X  *  265  nm  has 
melted,  the  second  harmonic  generation  shows  no  farther  increase.  The  results 
lead  to  an  energy  transfer  time  between  the  hot,  dense  carrier  plasma  and  the 
lattice  vibrations  of  less  than  two  picoseconds. 

A  careful  observation  and  analysis  of  the  photoelectric  emission  induced  by 
X  =  265  nm  picosecond  pulses  from  a  silicon  crystalline  surface  has  also  been 
completed.  The  data  show,  besides  the  ordinary  linear  photoelectric  effect,  an 
emission  proportional  to  the  square  of  the  intensity.  This  is  due  to  photoemis¬ 
sion  from  the  hot,  dense  carrier  plasma  created  by  the  same  pulse,  or  by  a 
preceding  picosecond  pulse.  The  lifetime  of  this  plasma  surface  layer  is  lim¬ 
ited  by  diffusion  of  the  hot  carriers  out  of  the  shallow  absorption  depth  for  X 
*  265  nm  radiation.  The  effect  of  space  charge  limitation  of  the  electron  emis¬ 
sion  in  an  intensity  regime  before  melting  occurs  has  also  been  demonstrated. 
After  the  onset  of  melting,  evaporation  of  atoms  and  ions  from  the  superheated 
fluid  silicon  takes  place.  The  positive  ions  reduce  the  space  charge,  and  the 
number  of  emitted  electrons  equals  that  of  the  emitted  positive  ions. 

Further  details  on  this  work  may  be  found  in  the  appendices,  which  are  pre¬ 
prints  or  reprints  of  technical  papers.  This  material  was  also  presented  at 
several  international  meetings.  Dr.  H.  Kurz  presented  invited  papers  at  the 
International  Conference  on  Quantum  Electronics  in  Anaheim,  California,  in  June 
1984,  and  will  also  present  an  invited  review  of  this  field  at  the  forthcoming 
meeting  of  the  Materials  Research  Society  in  Boston,  Massachusetts,  in  November 
1984. 


Recent  experiments  have  focused  on  the  problem  of  evaporation  of  neutral 


particles  from  the  superheated  fluid.  It  has  been  established  that  the  silicon 
readily  reaches  a  temperature  over  5000  E  at  a  fluence  of  500  mJ/cm  .  Thin  film 
SOS  samples  are  blown  off,  and  the  blow-off  phenomenon  has  hampered  the  study  of 
the  rate  of  evaporation. 

Theoretical  calculations  on  the  variation  of  the  effective  mass  of  the  car¬ 
riers  at  high  temperatures  and  densities  have  been  completed.  They  tend  to  con¬ 
firm  that  the  effective  mass  remains  much  smaller  than  the  free  electron  mass  up 
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to  densities  of  10  carriers/cc  and  up  to  the  melting  point. 

During  the  next  year  attention  will  be  focused  on  the  additional  damping  or 
energy  transfer  which  may  occur  when  the  plasma  frequency  of  the  carrier  gas 
exceeds  the  band  gap  energy.  It  should  be  possible  to  investigate  this  new 
regime  experimentally  in  germanium.  An  attempt  will  also  be  made  to  measure  the 
evaporation  rate  from  a  superheated  surface  layer  of  fluid  silicon  created  on  a 
bulk  single  crystal  of  silicon.  The  evaporated  atoms  will  be  collected  on  a 
transparent  sheet  and  their  density  measured  by  optical  means.  An  additional 
program  element  will  be  the  behavior  of  quantum  well  structures  at  high  carrier 
densities,  created  by  picosecond  pulses. 


B.  Picosecond  Pulsed  Later  Induced  Trent format  ions  in  Metals  and  Alloys 
The  personnel  associated  with  this  effort  were: 

F.  Spaepen,  Principal  Investigator 

W. K.  Wang,  Research  Fellow 

T.  Mizoguchi,  Visiting  Scientist  (summer) 

C.A.  MacDonald,  Graduate  Student 
S.M.  Prokes,  Graduate  Student  (eummer) 

L. C.  Chen,  Graduate  Student 

We  have  continued  the  program  on  the  formation  of  alloys  by  pulsed  laser 
heating,  mixing,  and  quenching.  We  have  concentrated  mainly  on  the  Nb-Si  sys¬ 
tem,  which  is  of  interest  for  its  superconducting  properties.  The  starting 
materials  were  made  by  ion  beam  sputtering,  in  our  recently  developed  apparatus 
(see  below,  and  attached  report),  and  consisted  of  three  types: 

1.  Homogeneous  Nb^Si^,  which  is  a  microcrystalline  Nb(Si)  solid  solution 
with  lattice  parameter  a  =  3.47  k. 

2.  Com  positionally  modulated  Nb^Si^j  (period  32  k) ,  consisting  of  fine 
grained  pure  Nb  and  amorphous  Si  layers. 

3.  Compositionally  modulated  (period  48  A),  of  the  same  structure 

as  2  above. 

After  laser  irradiation  (30  ps  TAG  or  30  ns  ruby  laser  pulses),  the  follow¬ 
ing  structural  changes  were  observed; 

1.  The  composition  modulation,  observed  by  x-ray  diffraction,  decayed  as  a 
result  of  large  fluence  irradiation. 

2.  The  Nb^gSijj  modulated  films  transformed  to  a  bcc  Nb(Si)  solid  solution 
with  lattice  parameter  a  *  3.30  k  after  30  ns  irradiation. 

3.  The  NbgjSi^p  modulated  films  transformed  to  an  as  yet  unidentified 


phase  (in  electron  diffraction)  after  30  ns  irradiation. 

4.  The  Nbgi Si^j  modulated  films  transformed  to  a  bcc  Nb(Si)  solid  solntion 
with  lattice  parameter  a  =  3.45  A  after  30  ps  irradiation. 

5.  The  Nh^Si^  films  transformed  into  a  mixture  of  an  a-15  type  Nb^Si, 
bcc  Nba  and  y-NbgSij  after  30  ns  irradiation. 

Some  preliminary  investigation  of  the  superconducting  properties  of  these 
materials  has  been  made.  It  has  become  clear  that  an  inductive  rather  than  a 
resistive  method  would  be  preferable  for  measurements  on  the  small  amounts  of 
transformed  material  available  after  laser  pulse  induced  transformation. 

Glass  formation  in  marginally  glass  forming  alloys,  which  formed  the  bulk 
of  last  year's  work,  was  continued  on  the  Fe-C  system,  in  which  so  far  no 
glasses  could  be  made  from  the  melt  by  conventional  quenching  methods.  Prelim¬ 
inary  investigations  show  that  30  ps  laser  quenching  leads  to  glass  formation  in 
alloys  containing  15  and  25  at.%  C. 

Ve  have  continued  our  study  on  the  direct  measurement  of  the  crystal  growth 
velocity  and  solute  trapping  in  pure  metals  and  simple  alloys.  The  pump-probe 
technique  has  been  improved  to  give  better  statistics.  We  have  also  begun  to 
investigate  interdiffusion  in  the  liquid  and  subsequent  solute  trapping  in  modu¬ 
lated  films,  using  Auger  spectroscopy  and  x-ray  diffraction. 

We  have  made  extensive  use  of  the  new  dual  gun  ion  beam  aputtering 
apparatus  developed  in  this  project  for  the  preparation  of  starting  materials. 

A  copy  of  a  technical  report  on  it  is  attached.  It  has  proved  very  versatile  in 
producing  both  homogeneous  alloys  and  various  compositionally  modulated  films. 
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Division  of  Applied  Sciences,  Harvard  University,  Cambridge,  MA  0213c,  USA 


Numerous  investigations  of  the  mechanisms  for  pulsed  laser-induced  phase  tran¬ 
sition  provide  ample  evidence  that  the  surface  of  metals  and  semiconductors 
undergoes  a  solid-liquid  phase  change  as  soon  as  a  critical  laser  fluence  is 
exceeded.  There  is  no  doubt  about  the  thermal  nature  of  the  phase  transition 
as  long  as  nanosecond  pulses  are  being  used  [1],  It  is  the  main  purpose  of 
this  contribution  to  clarify  whether  under  picosecond  irradiation  the  simple 
thermal  melting  approach  is  still  valid.  This  question  is  focused  on  the  time 
required  to  establish  thermal  equilibrium  between  carriers  and  phonons  [2], 

The  conventional  picosecond  pump  and  probe  technique  is  applied  [3,4].  Tne 
highly  focused  probe  beam  monitors  the  transmission  and  reflectivity  of  bulk 
surfaces  on  thin  silicon  films  on  sapphire  (SOS),  induced  by  an  exciting 
picosecond  pulse  at  0.532  pm.  The  wavelength  of  the  probe  pulse  is  varied 
between  0.532  and  2.8  pm.  By  comparing  the  changes  in  optical  properties  at 
different  wavelengths,  the  contributions  due  to  a  variation  in  lattice  tem¬ 
perature  TL  and  those  due  to  changes  in  carrier  density  N  are  separately 
determined  [5], 

At  the  doubled  frequency  of  a  Nd:YAG  laser  (hw  =  2.33  eV),  the  optical 
properties  of  picosecond-excited  silicon  are  mainly  determined  by  the  indirect 
gand  gap  transition.  The  phonon-assisted  indirect  absorption  a  depends 
strongly  on  the  lattice  temperature  Tl  [6],  In. SOS  samples  multiple  reflec¬ 
tions  from  air-silicon  and  silicon-sapphire  interfaces  enhance  the  optical 
detectivity  of  changes  Ac  in  the  real  part  of  the  dielectric  functions  con¬ 
siderably.  Time-resolved  analysis  of  the  reflectivity  signatures  of  SOS  sam¬ 
ples  reveal  the  interplay  between  thermal  (Ae  >  0)  and  free  carrier  (Ac  <  0) 
contributions  to  c  [9].  Due  to  Auger  recombination  the  plasma  density  is 
known  to  drop  below  1.2x10^  cm"3  after  200  ps.  At  this  time  the  plasma  con¬ 
tributions  become  negligible  at  0.532  ym  and  the  amount  of  lattice  heating 
can  be  determined  experimentally  by  comparing  samples  with  different  film 
thicknesses  [5]. 

In  Fig.  1  the  surface  temperature  developed  200  ps  after  the  exciting  pulse 
is  shown.  Picosecond  optical  probing  of  lattice  temperature  reveals  signifi¬ 
cantly  higher  values  than  those  derived  up  to  now  from  Raman  scattering  data 
using  ns  pulses.  Compared  to  bulk  surfaces  (F^h  =  200  mJ/cnr).  SOS  samples 
exhibit  a  lower  threshold  value  for  surface  melting  (160  mJ/cm^).  Because  of 
multiple  interferences  in  SOS  samples,  a  larger  amount  of  energy  densi  ty  is 
absorbed  in  the  silicon  film.  The  absorbed  energy  density  in  0.5  ym  thick 
silicon  films  is  determined  by  reflectivity  and  transmission  measurements  of 
the  pump  beam  (see  top  of  Fig.  1).  The  surface  temperature  increases  non- 

*Permanent  address:  C.E.N./Saclay ,  DPh.G/S.P.A. ,  91191  Gif-sur-Yvette  Cedex, 
France 
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Fig.  1.  Surface  temperatures  Ts 
versus  incident  laser  fluence  at 
10  and  200  ps  after  the  heating 
pulse  at  0.532  pm.  The  data 
are  evaluated  from  time-resolved 
reflectivity  and  transmission 
measurements  on  SOS  samples  using 
the  thermo-optic  data  of  ref.  [6] 


linearly  with  the  laser  fluence  F(2lj)  and  reaches  values  close  to  the  melting 
point  of  silicon  (Tm  =  1680  K)  at  a  fluence  level  (160  mJ/cnr)  where  the  re¬ 
flectivity  increases  rapidly  to  the  liquid  state  value.  Obviously  lattice 
heating  and  melting  occur  during  the  pump  pulse.  The  band  gap  shrinks  nearly 
instantaneously,  resulting  in  an  increase  of  the  pump  beam  absorbance.  This 
picture  is  confirmed  by  measuring  the  transmission  at  time  delays  shorter  than 
the  pulse  duration.  At  a  time  delay  of  At  =  10  ps,  significant  lattice  heating 
already  occurs.  The  result  shown  in  Fig.  1  is  obtained  by  temporal  deconvolu¬ 
tion  of  the  heating  and  probing  pulse.  Plasma  contribution  has  been  properly 
taken  into  account.  Clearly,  phonons  participating  in  the  indirect  absorption 
process  received  a  significant  amount  of  energy  deposited  primarily  in  the 
electron-hole  plasma.  The  observed  transmission  changes  are  in  excellent 
agreement  with  the  optical  heating  model,  where  an  instantaneous  energy  trans¬ 
fer  to  the  phonons  is  assumed. 

At  a  probing  frequency  below  the  indirect  band  gap,  free  carrier  contribu¬ 
tion  becomes  dominant.  Drude-like  plasma  resonances  have  been  observed  at 
2.8  urn  [7].  They  have  been  used  for  a  preliminary  estimate  of  plasma  densities 
far  below  the  threshold  for  melting.  At  this  probing  wavelength  the  plasma 
edge  is  reached  at  a  fluence  level  of  40  mJ/cm^  already.  Strong  free  carrier 
absorption  causes  transmission  drops  to  zero,  preventing  a  reliable  determina¬ 
tion  of  the  maximum  plasma  density  at  higher  fluence  levels.  For  this  pur¬ 
pose  thin  film  and  bulk  crystals  are  carefully  investigated  at  1.064  and  1.9  pm. 
At  these  probe  wavelengths  the  plasma  contributions  are  less  pronounced,  and 
measurements  up  to  the  melting  point  can  be  performed.  This  study  is  completed 
by  a  novel  three-pulse  technique,  where  the  first  pulse  at  0.532  ym  creates 
carriers,  a  second  pump  pulse  at  1.064  ym  couples  energy  to  the  plasma  by  free 
carrier  absorption  without  creating  new  carriers,  and  a  third  pulse  monitors 
the  induced  optical  changes  [8].  The  plasma  heating  of  the  second  pulse  does 
not  affect  the  susceptibility  associated  with  intraband  (Drude)  and  interband 
transitions.  Neither  the  reduced  mass  m*  nor  the  density  of  electron-hole 
pairs  is  changed  by  additional  heating  of  the  plasma.  The  probability  of  car¬ 
rier  multiplication  by  impact  ionization  is  very  low  [8,9].  The  energy  loss 
rate  of  the  carriers  to  the  phonons  is  faster  than  the  rate  for  impact  ioniza¬ 
tion,  which  is  expected  to  be  comparable  to  the  Auger  recombination  rate. 


As  a  selected  example,  the  transmission  of  bulk  silicon  at  1.064  ym  is 
plotted  versus  incident  laser  fluence  for  two  different  time  delays  At  in 


Fig.  2.  Transmission  of  bulk 
silicon  at  1.064  pm  versus  inci¬ 
dent  laser  fluence  at  0.532  pm 
for  two  different  time  delays 
between  pump  and  probe  pulse 


(J/tm*) 

Fig.  2.  At  At  =  200  ps  the  plasma  density  levels  off  to  a  constant  value 
above  F ( 2^, )  =  50  mJ/cm^  due  to  the  nonlinearity  of  the  Auger  recombination 
process.  The  lattice  temperature,  however,  increases  strongly  with  F(2oo). 
Under  these  conditions  the  free  carrier  absorption  at  1.064  pm  associated 
with  intraband  transitions  is  proportional  to  [m*  <t>  ( 2yt ) ,  where  m*  means 
the  reduced  optical  mass  for  electron-hole  pairs,  <t>  the  energy  averaged 
momentum  relaxation  time  of  the  carriers,  and  y  the  Auger  recombination  coef- 
fi  cient. 

The  solid  line  in  Fig.  2  represents  calculated  transmission  values  of  bulk 
silicon  using  y  =  2x10~3l,  m*  =  0.164  m0  and  averaged  momentum  relaxation 
time  <t>  =  7xl0-^5  x  300/Ti  ,  as  found  in  ns  experiments  [10].  The  agreement 
between  measured  and  calculated  data  is  excellent,  indicating  a  strong  phonon 
participation  in  the  carrier  momentum  relaxation,  even  at  plasma  density 
levels  at  which  electron-hole  scattering  dominates  the  relaxation  process. 

The  same  calculation  provides  satisfactory  results  for  the  transmission  at 
At  =  10  ps,  where  the  plasma  density  depends  on  the  laser  fluence  F(2w).  In 
this  case  the  plasma  density  has  been  calculated  with  the  standard  equations, 
including  ambipolar  diffusion  and  neglecting  impact  ionization.  The  calcula¬ 
ted  evolution  of  plasma  density  and  temperature  at  the  surface  of  silicon 
irradiated  with  increasing  laser  fluence  F(2uj)  are  finally  summarized  in 
Fig.  3.  These  calculated  data  are  cross  checked  with  experimental  data  at 
different  probing  wavelengths,  excitation  levels  and  time  delays.  They  have 
been  used  to  analyze  the  dependence  of  the  optical  reduced  mass  on  the  car¬ 
rier  density  in  reflectivity  measurements.  Excellent  agreement  with  the  data 
at  1.064,  1.9  and  2.8  pm  has  been  found  up  to  a  fluence  level  of  150  mJ/cm^, 
for  a  constant  optical  reduced  mass  of  (0.16  +  0.1 )m0.  The  stability  of  the 
Drude  term  indicates  a  "cold  plasma"  which  is  in  thermal  equilibrium  with  the 
lattice.  Close  to  the  threshold  value  for  melting,  the  maximum  plasma  density 
is  limited  to  1.4  xlO2'  cm’3,  mainly  due  to  Auger  recombination.  The  measured 
values  of  lattice  temperature  and  maximum  plasma  density  are  consistent  with 
optical  heating  and  thermal  melting  even  on  a  time  scale  of  picoseconds. 

This  research  was  supported  by  the  U.S.  Office  of  Naval  Research  under 
contract  N00014-83K-0030  and  by  the  Alexander  von  Humboldt  Foundation,  F.  R. 
Germany. 
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Fig.  3.  Calculated  plasma  densi¬ 
ties  and  temperatures  at  the  sur 
face  of  bulk  silicon  versus  inci 
dent  laser  fluence  at  0.532  pm 
using  the  optical  heating  model 
with  instantaneous  energy  trans¬ 
fer  during  the  picosecond  pulse. 
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ABSTRACT 

Three  different  photoelectric  regimes  are  obseived  in 
the  interaction  of  ii  ps,  266  nr.  laser  pulses  with  crystalline 
silicon  sam: les  versus  light  fluence.  A  superposition  of 
linear  and  quadratic  photoionizat ion  is  followed  by  a  space 
charge  limited  regime  u:  tc  the  critical  fluence  F^p  for  the 
surface  anorphization  where  highly  nonlinear  ion  emission  is 
observed.  Ior.  and  electron  emissions  become  equal  in  magnitude 
at  a  fluence  ~  2Fth*  T*le  aX3sence  of  observable  thermionic 
effects  indicates  that  thermal  equilibrium  of  the  electron- 
hole  plasma  and  the  lattice  is  reached  during  the  laser  pulse 
duration. 


INTRODUCTION 

The  central  question  in  the  investigation  of  pulsed  laser  annealing  of 
semiconductors  in  tne  picosecond  range  is  related  to  the  time-  scale  of  the 
energy  transfer  between  the  hot  electron-hole  plasma  created  by  the  laser 
pulse  and  the  crystal  lattice.  In  recent  times  much  experimental  effort  has 
been  devoted  to  the  study  of  this  energy  transfer.  This  goal  has  been  mostly 
pursued  through  optical  diaonostic  methods,  either  by  measuring  the  optical 
constants  of  the  senuconducotrs  with  rum*  and  probe  techniques  [1J  or  by 
direct  derivation  of  the  lattice  temperature  from  Stokes  and  anti-Stokes 
Ranan  scattering  (2]. 

Measurements  of  the  electron  and  ion  emission  from  the  semiconductor 
surface  during  the  interaction  with  picosecond  laser  pulses  offers  an  alter¬ 
native  approach  to  the  problem.  Here,  one  probes  the  material  only  in  its 
outermost  layers,  and  therefore  the  data  are  insensitive  to  the  gradients  in 
the  profiles  of  plasma  density  or  lattice  temperature.  Thermionic  emission 
is  extremely  sensitive  to  the  electron  temperature.  Conversely,  the  emission 
of  positive  ions  is  an  indicator  of  a  hot  lattice. 

In  this  paper  detailed  measurements  of  the  charges  emitted  by  crystal¬ 
line  silicon  samples  irradiated  with  266  nm,  15  ps  laser  pulses  are  reported. 
This  work  extends  the  measurements  of  Liu  et  al.  (3]  to  laser  fluences  about 
four  orders  of  magnitude  below  the  critical  value  for  a  surface  phase  trans¬ 
formation.  One  of  the  main  results  of  this  analysis  is  the  absence  of  ob¬ 
servable  thermionic  effects  throughout  the  whole  range  of  laser  fluence. 

Thus,  an  upper  limit  to  the  average  electron  temperature  of  «  3000  K  within 
the  laser  pulse  can  be  deduced. 

EXPERIMENTAL  HETHCD 

The  source  of  excitation  of  the  silicon  samples  is  a  15  ps,  266  nm  laser 
pulse  obtained  by  frequency  quadrupling  the  output  of  a  Nd:YAG  30  ps  laser 
pulse.  The  radiation  is  focused  by  quartz  optics  into  a  vacuum  chamber  in 
which  the  silicon  sample  is  located.  The  angle  of  incidence  is  48° ,  and  the 
reflected  portion  of  the  beam  escapes  from  the  chamber  through  a  second 
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window.  A  1  r  diameter  wire  placed  at  2.5*3  mm  from  the  sample  surface- 
acts  as  a  charge  collector.  Botn  sample  and  wire  are  electrically  msulatvd 
from  the  body  of  the  vacuum  chamber  and  can  be  independently  biased  up  tr 
«*•  4  kV.  Two  configurations  for  the  measurements  are  possible.  One  can  mea¬ 
sure  either  the  charges  collected  on  the  wire  or  the  ones  escaping  from  the 
sample.  Both  methods  have  been  shown  to  give  similar  results.  Most  of  the- 
measurements,  however,  have  been  performed  by  measuring  the  total  charges 
escaping  from  the  target.  A  charge  amplifier  connected  to  the  sample  col¬ 
lects  the  charges  and  provides  a  signal  proportional  to  the  collected 
charges  to  the  data  acquisition  system.  The  calibration  of  the  amplifier 
has  been  checked  for  high  charge  signals  against  the  charge  signal  obtained 
at  the  oscilloscope  through  a  coaxial  cable  of  known  capacitance.  The  read¬ 
ings  agree  within  a  factor  1.5.  An  absolutely  calibrated  photodiode  and 
neutral  density  filters  provided  the  total  energy  value  of  each  laser  pulse. 
Laser  and  charge  data  were  directly  plotted  by  an  automatic  data  acquisition 
system. 

Measurements  of  the  spatial  beam  profile  on  the  target  were  performed 
for  each  experimental  run  by  evaluating  the  area  of  the  amorphous  spots  at 
various  laser  fluences,  as  described  in  reference  [4] . 

The  experiments  were  performed  in  a  vacuum  of  5  * 10~7-5  *  10"8  Torr. 

The  samples  were  Si  (111)  and  Si  (100)  lightly  doped  p-samples  with  1-5  fi-cm 
and  15-50  respectively . 

RESULTS  AND  DISCUSSION 

The  electron  and  ion  emission  from  Si  (111)  and  Si  (100)  has  been  mea¬ 
sured  in  many  experimental  runs  with  different  focusing  conditions.  The 
(1/e)  radius  p  of  the  laser  intensity  profile  on  the  sample  has  been  varied 
between  12  and  178  Um.  The  electron  signal  does  not  depend  on  the  contami¬ 
nation  of  the  sample  surface.  Cleaning  procedures  r.uch  as  multiple  shot  ex¬ 
posure  below  the  fluence  for  morphological  surface  damaqe  (5)  do  not  modify 
the  observed  features  in  the  collected  charges.  Contrary  to  expectation,  no 
qualitative  difference  in  the  emission  has  been  found  between  Si  (111)  and 
Si  (100)  oriented  samples,  despite  the  differences  in  ionization  potentials 
(6,7).  Presumably  the  surface  oxide  layer  present  on  the  samples  exposed  to 
air  is  not  removed  by  laser  pulse  cleaning  and  thus  lowers  the  ionization 
potential  at  the  surface  of  the  sample  (8) . 

Figure  1  shows  the  total  of  collected  electron  and  ion  charges  versus 
laser  fluence  in  a  typical  plot  for  Si  (111).  The  behavior  of  the  electron 
emission  can  be  divided  into  three  regimes.  At  lower  fluences  (region  I)  we 
have  a  progressively  nonlinear  departure  from  a  pure  linear  (slope  =  1)  re¬ 
lationship  between  electron  signal  and  laser  fluence  This  indicates 

that  the  energetically  allowed  single  photon  effect  is  accompanied  by  a  non¬ 
linear  photoelectric  effect.  In  a  second  range  of  laser  fluences  (region  II) 
a  decrease  in  the  slope  of  the  curve  is  observed.  The  collected  charge  sig¬ 
nal  is  here  dependent  upon  the  bias  voltage  and  scales  approximately  with 
Vbias*  ^is  behavior  is  then  interrupted  at  a  fluence  Fjjj  =  0.08  J/cm2, 
where  amorphous  spot  formation  and  positive  ion  emission  occur  simultane¬ 
ously.  Above  F^jlJ  (region  III)  a  steep  increase  of  the  electron  emission 
occurs.  In  this  regime  each  experimental  point  is  obtained  on  a  fresh  por¬ 
tion  of  the  sample.  Ion  emission,  on  the  other  hand,  is  characterized  by  a 
very  steep  behavior.  There  is  no  significant  amount  of  ions  emitted  in  re¬ 
gions  I  and  II.  The  ion  signal  emerges  from  the  noise  at  about  one  quarter 
of  the  critical  fluence  F^.  The  first  part  of  the  ion  emission  has  a  rela¬ 
tively  low  slope.  It  is  followed  by  a  jump-like  behavior.  Within  25%  of 
the  critical  value  fJJJ  the  emission  increases  by  more  than  two  orders  of 
magnitude.  At  the  critical  fluence  for  amorphous  formation,  the  curve 

bends  to  a  lower  slope,  and  equal  amounts  of  electrons  and  ions  are  observed 
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FIG.  1.  Era. t ted  cnarae  versus- 
incident  laser  fluer.ce  for  (111; 
oriented  silicon  samples.  To*, 
laser  spot  p  at  1/e  is  5"  Ur.. 
Full  dots,  electrons;  triangles, 
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above  in  region  III. 

A  more  illustrative  way  of  displaying  these  three  regimes  of  electron 
emission  is  shown  in  Fig.  2 ,  where  the  ratio  Q/ Zn  of  the  total  electron 
charge  Q  to  the  incident  laser  energy  Er  (in  Joules)  is  plotted  versus  laser 
fluence  F(J/cm*-).  The  behavior  at  low  fluences  (region  I)  up  to  laser 
fluences  ~  F^/15  is  linear,  as  shown  on  an  enlarged  scale  in  Fig.  3.  This 
indicates  that  a  superposition  of  linear  and  quadratic  ionization  processes 
is  effective  m  this  range  of  laser  fluences.  The  functional  dependence  of 
the  current  density  J(r,t)  versus  laser  intensity  I(?,t)  is  described  by 


where  e  is  the  electron  charge,  hv  the  photon  energy  and  K^,  1*2  t^e  co effi¬ 
cients  for  linear  and  quadratic  photoemission ,  respectively.  Integrating 
this  equation  over  the  space  and  time  with  the  Gaussian  beam  profile  results 


n  2*2  (hv)  \ 

Here  T  is  the  laser  pulse  duration,  £n  the  total  energy  of  the  laser  pulse 
(J)  and  F  its  fluence  (J/cm2). 

The  value  of  the  quantum  yield  for  single  photon  emission  at  4.66  eV , 
as  deduced  from  the  experiment,  equals  8 (*  2)  * 10“7  electrons  per  absorbed 
photon,  in  agreement  with  data  found  in  the  literature  for  silicon  samples 
exposed  to  air  (7J .  The  experimental  value  for  the  quadratic  effect,  3  + 

1.5  * 10“32  cm^  s  is  consistent  with  the  results  obtained  by  Bensoussan  et  al . 
[9) ,  if  extrapolated  to  higher  photon  energies.  This  nonlinear  process  com¬ 
pletely  dominates  the  electron  emission  over  a  wide  range  of  laser  fluences. 
This  can  be  attributed  to  a  much  higher  escape  probability  for  the  two-photon 
excited  electrons  which  have  an  intrinsically  higher  (up  to  ~  4.6  eV)  kinetic 
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fig.  2.  The  data  are  replotted  ir.  terns  FIG.  3.  Ratio  Q/Zr  of  emitted 

of  the  ratio  Q/Z n  (emitted  negative  charge  negative  charge  to  incident 

over  incident  laser  energy)  versus  inci-  laser  energy  versus  incident 

dent  laser  fluence  for  Si  (111).  laser  energy  at  different  focus¬ 

ing  conditions  on  a  Si  (ill) 
sample.  The  values  of  the 
laser  spot  radii  are  shown  for 
each  curve . 

energy  than  the  electrons  emitted  by  a  single  photon  process,  which  have 
about  0.06  eV  kinetic  energy. 

The  nonlinear  behavior  at  fluences  higher  than  ~  f£k/15  is  masked  by 
space  charge  limitations.  This  regime  of  emission  depends  on  the  bias  volt¬ 
age  on  the  collecting  electrode  and  cannot  be  neutralized  at  the  highest  bias 
voltages  permissible  in  our  experimental  conditions.  Space  charges  rapidly 
neutralize  the  applied  fields  at  the  surface.  The  electron  emission  is 
limited  to  the  leading  edge  of  the  picosecond  pulse. 

Simultaneously  with  the  observation  of  amorphous  spot  formation  at  the 
critical  laser  fluence  F^£,  sharp  rises  in  ion  and  electron  emission  occur. 

The  ion  emission  has  a  step-like  behavior,  and  the  electron  space  charge 
limitation  is  suddenly  released.  These  discontinuities  show  that  the  surface 
conditions  of  the  sample  are  suddenly  modified  at  this  fluence,  as  one  would 
expect  for  a  phase  transition  to  the  molten  state  of  silicon. 

Thermionic  emission  is  a  very  steep  function  of  the  electron  temperature 
Te .  Therefore,  it  should  manifest  itself  with  a  highly  nonlinear  behavior 
even  in  region  I.  In  Fig.  3  the  ratio  Q/En  versus  incident  laser  energy  En 
for  various  laser  spot  sizes  in  region  I  shows  that  the  behavior  is  strictly 
linear.  Thus  the  thermionic  contribution  below  the  saturation  equals  at  most 
the  accuracy  of  the  data.  Application  of  the  Richardson-Dushman  equation 
gives  an  upper  limit  of  ~  2200  X  for  the  average  electron  temperature  during 
the  laser  pulse  in  region  I.  In  region  III  the  equal  amounts  of  electrons 
and  ions  observed  above  2F*JJ  reduce  the  possible  thermionic  emission  contri¬ 
bution  to  less  than  10“*2  Coulombs  in  this  regime.  This  gives  an  upper  limit 
of  ~  3000  K  for  the  average  electron  temperature  during  the  phase  transition. 

CONCLUSIONS 

The  results  of  our  measurements  can  be  summarized  as  follows.  The  elec¬ 
tron  charges  collected  at  lower  laser  fluences  must  be  attributed  completely 
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to  trie  sut-cr:  rsitic:.  of  single-photon  photoemission  and  a  quadratic  phot  o¬ 
electric  effect.  No  ex;.K. nner.tal  evidence  for  thermionic  emission  is  found 
in,  our  measurements.  _T:.is  provides  an  upper  limit  of  303o  K  to  the  average 
electron  temperature  Tfe  during  the  laser  pulse  duration.  The  data  are  con¬ 
sistent  with  the  establishment  of  equilibrium  between  the  carrier  and  the 
lattice  temperature  during  the  laser  pulse  {15  ps).  A  space  charge  limited 
regime  is  observed  for  the  electron  emission  at  fluences  between  F^'15  and 
Sufficient  charge  densities  (Coulonbs/cn^)  can  build  u:  to  screen  the 
applied  field  lr.  a  fraction  of  the  laser  pulse  duration.  At  the  critical 
fluence  for  the  phase  transformation ,  melting  resulting  in  subsequent  forma¬ 
tion  of  an  amorphous  pattern,  the  sudden  release  of  the  space  charge  regime 
and  the  appearance  of  positive  charges  are  consistent  with  the  notion  of 
surface  melting  during  the  laser  pulse  110]. 
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ABSTRACT 

Time-resolved  studies  cf  reflectivity  and  transmission 
at  C.532  ur.,  1.0c4  ».m  and  2.t>  ..t  of  thin  silicon  films  fol¬ 
lowing  irradiation  with  ps  pulses  at  0.532  urn  have  been  per¬ 
formed.  The  formation  cf  the  electron-hole  plasma  and  the 
evolution  of  lattice  temperature  is  investigated  as  a  function 
of  pum:  fluer.ce  and  time  delay.  <2uantitative  determination 
of  the  plasma  densities  and  lattice  temperature  up  to  the 
melting  temperature  shows  that  the  maximum  plasma  density  is 
limited  to  ~  1  w  10*A  cm“3  by  Auger  recombination  even  on  a 
time  scale  of  picoseconds  at  fluences  sufficient  to  cause  th»_- 
phas*  transition.  The  thermal  nature  of  the  phase  transition 
is  confirmed. 


Numerous  investigations  of  the  change  m  optical  properties  of  silicon 
have  beer,  carried  out  as  a  function  of  picosecond  irradiation  fluence  [1-4;  . 
These  data  indicate  melting  of  the  surface  when  the  pulse  fluence  exceeds  a 
critical  threshold  value.  This  paper  shows  that  more  quantitative  informa¬ 
tion  about  the  lattice  temperature  and  the  electron-hole  plasma  kinetics  can 
be  obtained  by  varying  the  wavelength  of  the  probing  beam  and  its  time  delay 
with  respect  tc  the  heating  pulse. 

The  real  part  of  the  dielectric  function  at  a  probing  frequency  u.  below 
the  direct  bandqa; , 

C'  =  n^(T,N)U  -  k2(T,N>]  -  AN  a: 


where 


A  = 


w2  W  Bih  / 


increases  with  the  lattice  temperature  T  due  to  indirect  mterband  transi¬ 
tions  and  decreases  with  the  number  of  electron-hole  pairs  N.  The  optical 
reduced  mass,  m*  =  (m~*  ♦  m^)~*  may  depend  on  the  structure  in  the  vicinity 
of  the  conduction  band  minimum  and  valence  band  maximum  covered  by  the  ener¬ 
getic  distribution  of  the  carriers  during  the  time  of  observation. 


The  imaginary  part  of  the  dielectric  function 


£"  -  2n^(T,N)kL<T,N)  +  BN 

where 


(2 


B  = 


4 ire2  (  1 

3  Im  'T  > 
uj  \  e  e 


VV 


) 


increases  with  the  temperature  T  and  plasma  density  N. 
parameter,  the  scattering  times  of  the  carriers  (Te  and 
their  energy  distribution  come  into  play. 


As  a  third  unknown 
Ih)  averaged  over 


At  a  wavelength  of  532  nm  the  real  part  e*  reveals  the  opposing  effects 
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of  increasing  lattice  ten:  eraturv  and  increasing  ;  lasna  density.  Trie  ina  .Ji¬ 
nary  fart  der-er.ds  mainiv  or.  the  lattice  temperature ,  while  plasma  abscr;  ti 
is  negligible  at  5.31  nr-.  (5,61.  By  contrast,  at  2 .  b  u~  and  1.9  Ur  both  t*  and 
c“  are-  dominated  by  plasma  effects  (7). 

Careful  comparison  between  the  time-resolved  spectra  of  the  reflectivity 
and  transmission  at  these  wavelengths  permits  the  quantitative  determinaticr. 
of  the  lattice  temperature  and  plasma  density  during  and  after  ps  excitation*, 
by  photons  whose  energies  are  about  twice  as  large  as  the  indirect  band  gaj 
of  Si , 

To  avoid  undesirable  diffusion  phenomena  and  obtain  uniform  temperature 
and  t lasna  profiles  within  the  penetration  depth  of  the  exciting  i urn;  laser, 
the  measurements  have  been  performed  or.  thin  (C.5  urn)  films  of  silicon  on 
sapphire  (SOS).  In  addition,  the  multiple  interferences  between  air-silicon 
and  silicon-sapphire  interfaces  increase  the  sensitivity  to  changes  of  the 
dielectric  function. 

Further  understanding  of  the  picosecond  dynamics  of  semiconductor  car¬ 
riers  and  lattice  has  been  obtained  by  a  novel  three-pulse  technique.  The 
first  pulse  at  532  nm  creates  carriers.  A  second  pump  pulse  at  1.004  oi 
30  ps  duration  adds  energy  to  the  plasma  by  free  carrier  absorption,  without 
changing  its  density.  The  effect  of  this  additional  plasma  heatinc  is 
studied  by  a  third  probing  pulse  at  532  nr,  at  different  time  delays  with  re¬ 
spect  to  the  first  two  heating  pulses  [8] .  If  the  excess  energy  is  retained 
m  the  plasma  for  times  comparable  to  Auger  recombination  time,  the  tempera¬ 
ture  of  the  electron-hole  plasma  is  increased  above  the  band  gaj  .  Impact 
ionization  and,  consequently,  a  lowering  of  the  dielectric  function  are 
expected . 

Tne  experimental  set-up  for  all  these  measurements  is  similar  to  that 
described  in  reference-  (3|.  Due  to  Auger  recombination  the  electron-hole 
plasma  density  drops  below  1  *  1020  cm-3  after  150  ps ,  and  the  plasma  contri¬ 
bution  to  £'  and  C"  becomes  negligible  for  532  nm  liqht.  At  a  time  delay  of 
200  ps  the  lattice  temperature  can  be  directly  measured  by  transmission 
changes  of  silicon  film  at  various  fluence  levels.  The  SOS  samj le  is  heated 
by  a  20  ps  pulse  with  a  diameter  of  3C0  *jr..  The  probe  pulse  at  the  same 
wavelength  is  focused  to  30  urn  and  aligned  directly  at  the  center  of  the  pun: 
spot.  The  polarization  of  the  two  beams  is  orthogonal  to  suppress  undesir¬ 
able  interferences.  The  fluence  of  the  probe  pulse  is  kept  at  0.1%  of  the 
pump  fluence  to  avoid  reheating  phenomena. 

As  shown  in  Fig.  1,  the  SOS  film  exhibits  significant  changes  of  reflec¬ 
tivity  and  transmission  with  increasing  fluence  level  of  the  pum$  pulse  F^. 
The  reflectivity  increases  gradually  from  15%  to  40%,  followed  by  a  flat 
level  up  to  0.16  J/cm* .  Above  this  level  a  sharp  rise  to  69%  occurs,  which 
corresponds  to  the  metallic  reflectivity  of  liquid  silicon  at  an  anqle  of  26° 
as  used  in  this  experiment.  The  flat  level  portion  reflects  the  conditions 
of  bulk  reflectivity  because  of  the  nearly  complete  suppression  of  the  mul¬ 
tiple  interferences  by  the  strong  increase  of  the  film  absorbance. 


The  transmission  drops  continuously  from  an  initial  value  of  48%  to  the 
detection  limit  of  1%  at  0.2  J/cm2.  This  limit  corresponds  to  an  average 
optical  absorption  a  *  4^nk/X  >  7*  104  cm-*,  nearly  one  order  of  magnitude 
above  the  room  temperature  value  of  a  ~  104  cm-*.  According  to  the  measure¬ 
ments  of  Jellison  and  Modine,  the  absorption  coefficient  a  increases  exponen¬ 
tially  with  temperatures  a  *  a0exp(T/6),  where  8  =  430  *  30  K  and  aQ  is  a 
constant  factor  19 J .  We  determined  aQ  *  4200  cm"*  for  the  GOG  sample  at  room 
temperature.  From  the  transmission  data  of  Fig.  1  we  extract  the  surface  and 
average  temperatures  in  the  following  way.  We  calculate  the  temperature  pro¬ 
file  T  (x )  after  200  ps  for  various  laser  fluences,  taking  into  account  the 
temperature  dependence  of  optical  absorption,  thermal  conductivity  and  heat 


capacity.  These  T(x)  profiles  are  converted  into  absorption  profiles  a(x)  = 
aoexp(T(x)/0J  and  compared  with  the  optical  densities  o  *  d-*  a(x)dx 
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FIG.  1.  Reflectivity  and  transmis¬ 
sion  at  y  =  53.i  nr.,  probed  with  a 
delay  At  =  201  ps  after  the  pump 
pulse  at  y  -  5j 2  nm,  versus  incident 
fluence  in  a  0.5  'urn  thick  SOS  sample. 


obtained  from  the  transmission  data. 

The  corresponding  surface  tempera¬ 
tures  Ts  and  average  temperature  T  = 
d~A  T(x)dx  of  SOS  are  plotted  in 
Fig.  2  versus  the  incident  laser 
fluence  Fq1".  It  should  be  noted  that 
the  multi; lo  interferences  of  SOS 
lead  to  rapidly  varying  spatial  modu¬ 
lation  of  tne  energy  absorption  during 

the  ps  irradiation.  Therefore,  we  measured  in  addition  the  reflectivity  R 
and  transmission  T  at  a  time  delay  At  =  0.  From  the  fraction  of  absorbed 
energy  A  =  1  -  R  -  T,  we  calculated  the  expected  temperature  rise  using  the 
simple  thermal  model.  There  is  excellent  agreement  between  the  two 
approaches . 


FIG.  2.  Average  temperature  T  and 
surface  temperature  Ts  in  SOS  sample, 
obtained  from  the  transmission  data 
versus  incident  fluence. 


This  optical  probing  of  the  lattice  temperature  reveals  significantly 
higher  values  than  those  derived  from  Raman  scattering  (10]  .  When  our  data 
are  extrapolated  to  higher  fluence  levels  along  the  dashed  lines  m  Fig.  2, 
the  melting  temperature  Of  1680  K  is  reached  at  the  surface  of  SOS  for  a 
laser  fluence  above  160  mj/cm2.  It  is  precisely  at  this  fluence  level  that 
the  reflectivity  of  SOS  starts  to  rise  to  the  liquid  state  value,  as  shown  in 
Fig.  1. 

At  shorter  time  delays  At  between  pump  and  probe  pulses,  the  negative 
plasma  contribution  (“AN  =  -  (4"Te Vcu*m#)N)  competes  with  the  positive  contri¬ 
bution  from  the  increased  lattice  temperature.  The  complete  time-resolved 
spectra  of  reflectivity  and  transmission  of  SOS  shown  in  Fig.  3  reveals  the 
temporal  change  of  the  dielectric  constants  6 ’  and  e’*  by  the  formation  of 
electron-hole  plasma  and  evolution  of  the  lattice  temperature.  The  trans¬ 
mission  drop  is  obviously  terminated  immediately  after  the  excitation  pulse, 
indicating  that  the  lattice  heating  is  completed  at  this  time  and  the  phonons 
contributing  to  the  indirect  absorption  process  have  reached  a  steady  state 
temperature . 

The  reflectivity,  however,  rises  to  a  maximum  at  At  —  0  ps  and  passes 
through  a  minimum  where  the  plasma  contribution  AN  is  compensated  by  the  tem¬ 
perature-induced  increase  of  n^fT).  The  slow  recovery  of  the  reflectivity  to 
the  final  value  of  22%  is  clearly  due  to  decay  of  the  plasma  density  via 
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Flo.  3.  Time -re solved  reflectivity 
arid  transmission  in  SOS  induced  by  , 
l  un:  pulse  at  53 2  nm  with  52  cr*" 
and  probed  at  the  same  wavelength  a‘ 
different  time  delays  (solid  line'. 
Trie-  dotted  line  indicates  changes 
induced  by  a  second  pump  bear  at 
1004  nm  with  25  mJ/cm^  applied  with 
the  time  delay  At =  100  ps  after 
the  first  pump  pulse  at  532  nr. 

(50  mj/cm2 ) . 


Auger  recombination.  To  calculate  the  oiasna  density,  *  role  of  the-  im;  act 
lu.Mzatior.  during  the  laser  irradiation  has  to  be  clanf.i-.-d. 

A  first  pulse  at  532  nm  with  50  mJ/cm"  creates  a  carrier  plasma  and  ele¬ 
vates  the  lattice  tenj>erature  (full  line  in  Fig.  3).  A  second  pulse  with 
30  p?  duration  at  1.064  un  adds  around  1  eV  excess  energy  to  the  plasma  at  a 
time  delay  =  100  ps.  A  third  pulse  at  532  nr.  probes  the- resulting 

changes  in  reflectivity  AR  and  transmission  AT  at  a  time  delay  At 31.  The 
dotted  line  in  Fig.  3  shows  the  effect  of  the  second  heating  pulse  (25 
mJ/cn2).  A  further  increase  of  lattice  temperature  is  observed,  but  no  evi¬ 
dence  for  a  change  in  the  plasma  term  4TTe^N/u»^m*  can  be  deduced.  Obviously 
neither  the-  reduced  mass  m*  nor  the  density  of  electron-hole  pairs  is  changed 
by  the  additional  heating  of  the  plasma.  As  in  the  case  of  532  nm  excita¬ 
tion,  the  optical  input  energy  is  immediately  transferred  to  the  phonon 
system.  The  plasma  cools  off  so  rapidly  that  the  probability  for  the  crea¬ 
tion  of  new  electron-hole  pairs  is  very  low. 

Thus  the  electron-hole  plasma  kinetics  is  mainly  governed  by  the  optical 
generation  rate  G,  the  Auger  recombination  rate  \N3  and  the  diffusion  rate 
»N0a  where  Da  is  the  ambipolar  diffusion  coefficient. 

The  three-pulse  plasma  heating  technique  permits  a  precise  determination 
of  the  temporal  dependence  of  the  free  carrier  absorption  &pCA  =  nofca'  by 
measuring  the  additional  amount  of  lattice  heating  at  different  delay  time- 
In  Fig.  4  the  probing  pulse  is  set  at  a  time  delay  of  At  =  400  ps , 
monitoring  the  relative  change  of  reflectivity  AR  and  transmission  AT. 

Clearly  the  additional  heating  of  the  lattice  is  proportional  to  the  number 
of  electron-hole  pairs  N.  The  heating  peaks  at  a  time  delay  At21  of  10  ps , 
indicating  the  moment  when  the  plasma  density  reaches  a  maximum. 

The  decrease  of  AR  and  AT  at  larger  At2^  reflects  again  the  decay  of  the- 
plasma  density  via  Auger  recombmat ion.  This  thermometric  determination  of 
plasma  relaxation  is  by  far  more  sensitive  than  direct  optical  measurements. 
The  solid  line  in  Fig.  4  shows  the  calculated  relative  reflectivity  and 


FIG.  4.  Relative  changes  of  reflec¬ 
tivity  (AR)  and  transmission  (AT) 
probed  at  a  delay  fixed  at  400  ps 
induced  by  a  second  pulse  at  1064  nm 
with  a  fluence  of  40  mJ/cm2  as  a 
function  of  its  time  delay  At2^. 


transmission  changes  cf  The  plasma  density  N;x,t)  has  been  calculated 

numencallv  with  >  =  4  '  l?-"'*  cr^/sec,  D,  =  6  •  10 3  T7*  crt/s  (20  '  3C  ’T , '» 

T  “  L>  L 

and  Cj. =  5  *  1J  7 .  '32.  cx"  .11].  The  calculation  assumes  that  the  ab¬ 
sorbed  laser  energy  is  distributed  immediately  to  the  lattice.  The  relative- 
changes  of  R  and  T  are  evaluated  by  using  conventional  thin  film  optics  ex¬ 
pressions  at  a  fluence  level  of  F--*-  =  50  mJ  ''em- .  The  plasma  density  reaches 
a  maximum  of  ~  4  *  102*  cm-3  at  2>t  »  10  ps  and  drops  to  6.5  *  10*'  cm''  at 
300  ps,  exactly  reproducing  the  relaxation  by  Auqe-r  recombination. 


By  changing  the  probe  wavelength  to  2.8  um,  the  contribution  of  the  lat¬ 
tice  temt'erature  to  C’  and  £"  becomes  negligible  compared  to  the  plasma  ef¬ 
fects.  The  transients  in  reflectivity  and  transmission  of  SOS  are  exclu- 
dively  dominated  by  the  plasma  terms  AN  and  BN  at  this  wavelength.  Probed  at 
a  time  delay  of  1C  ps,  where  the  plasma  density  reaches  a  maximum,  the  trans¬ 
mission  of  SOS  drops  to  zero  for  2.6  urn  already  at  low  laser  fluences  of  the 
excitation  beam,  as  shown  in  Fig.  5.  The  initial  rise  in  the  transmission 
below  20  mj''cm“  is  solely  due  to  interference  effects  in  SOS  samples. 


The  reflectivity 
density  exceeds  Nc  =  r 


exhibits  a  plasma  resonance,  indicating  that  the  plasma 
(c/eX)*  =  2.5*  ir:23  — ~3 
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where  C-, 


11.6  is  the 


optical  reduced  susceptibility  masses  of  electron-hole  plasma  in  equilibrium 
with  the  lattice  phonons.  As  soon  as  the  plasma  edge  is  reached,  the  trans¬ 
mission  drops  to  zero  and  no  further  probing  of  the  plasma  kinetics  is 
possible . 


To  determine  the  plasma  densities  at  higher  fluence  levels  close  to  the 
threshold  for  the  phase  transition,  probing  at  1.064  Um  is  preferable.  At 
this  wavelength  the  reflectivity  does  not  exhibit  a  Drude-like  resonance. 

The  density  appears  to  be  limited  below  the  critical  density  Nc  =  1.76  *  lC*"4- 
cm”3  for  this  wavelength.  The  transmission  remains  measurable  up  to  the 
threshold  fluence  for  melting.  Again,  the  slight  increase  of  the  transmis¬ 
sion  at  low  fluences  is  due  to  multiple  interferences  of  the  probing  light 

between  the  SOS  interfaces.  Around  F2~  =  0.160  J/cm2,  where  the  SOS  sample 

o 

starts  to  melt,  transmission  drops  to  34i.  With  the  free  carrier  absorption 
O  =  5.1  *  10”*8  T/3oj  cm-  reported  in  the  literature  and  with  our  numeri¬ 
cally  calculated  distribution  functions  for  N(x)  and  T(x) ,  a  density  at  the 
surface  Kfi  =*  1  *  102*  cr.“3  is  derived.  This  value  is  in  good  agreement  with 
theoretical  calculations  of  the  spatial  and  temporal  development  of  the 
plasma  density  N(x,t)  based  on  thermal  equilibrium  between  plasma  and  lattice 
phonons . 


In  conclusion,  time-resolved  measurements  of  reflectivity  and  transmis¬ 
sion  of  SOS  samples  as  a  function  of  pumf  fluence,  time  delay  and  of  probing 
wavelength  provide  quantitative  data  about  the  lattice  temperatures  and 
plasma  density  during  and  after  the  irradiation  with  picosecond  pulses. 

Close  to  fluence  threshold  for  phase  transitions  at  the  sample  surface,  the 
lattice  temperature  reaches  the  melting  point  and  the  plasma  density  peaks  to 


Fig.  5.  Time-resolved  transmission 
in  SOS  (0.5  Um)  induced  by  a  pump 
pulse  at  532  nm  at  different 
fluences  and  probed  at  1.064  and 
2.8  um  at  a  time  delay  of  10  ps . 
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12^ cr~~ ,  limited  r-y  Auger  reconi.ir.at.ior. .  There  is  nc  evidence  for  ir;  a:: 
ionization  throughout  the  range  cf  laser  f luc-r.ces ,  indicating  a  fast  energy 
loss  of  the  electron-hole  plasma  via  phonon  emission.  Tne  characteristi c 
energy  transfer  tine  is  less  than  a  few  picoseconds,  m  agreement  with  tr.t 
theoretical  expectations  based  on  deformation  potential  scattering.  The  lat¬ 
tice  reaches  the  melting  point  during  the  picosecond  pulse,  if  the  energy 
fluence  exceeds  the  threshold  value  necessary  to  observe  a  phase  transition. 

This  research  was  supported  by  the  United  States  Office  of  Naval 
Research  under  contract  N0G014 -63-K-0030 . 
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Abstract.  Picosecond  time-resolved  reflectivity  measurements  on  bare  silicon  surfaces  and 
silicon  surfaces  with  oxide  layers  reveal  very  fast  heat  diffusion  and  material  evaporation  on 
subnanosecond  time  scales.  With  a  thick  oxide  layer  resolidification  of  a  molten  silicon 
surface  can  take  place  in  a  few  hundred  picoseconds.  At  high  laser  fluences,  vaporization 
processes  take  only  a  couple  of  lOOps. 

PACS:  64.70.Fx.  44,  79.20.Ds 


Recent  experimental  results  [1,  2]  on  ultrafast  pulsed 
laser  interactions  with  silicon  have  demonstrated  that 
thermal  equilibrium  between  the  dense  photo-excited 
electron-hole  plasma  and  the  silicon  lattice  is  es¬ 
tablished  on  a  picosecond  time  scale.  Experiments  with 
femtosecond  laser  pulses  [3]  indicated  subpicosecond 
energy  transfer  times  from  the  carriers  to  the  lattice 
and  subsequent  structural  changes  of  the  silicon  lattice 
in  a  few  picoseconds.  Our  earlier  experiments  [4] 
measured  on  picosecond  time  scales  the  temperature 
rise  of  a  silicon  lattice  induced  by  20-ps  laser  pulses. 
When  the  laser  pulse  has  sufficient  fluence  to  induce 
melting,  heating  rates  can  be  as  high  as  1014  C/s. 
Transient  lattice  heating,  melting,  and  overheating  of 
the  melt  are  observed  during  the  20-ps  pulse  duration 
[5].  The  extremely  high  heating  rates  on  the  surface 
create  very  steep  gradients  in  the  temperature  profiles 
which  then  induce  very  high  cooling  rates  through  fast 
diffusion  of  the  heat  into  the  bulk  substrate  after  the 
pulse.  These  ultrafast  heating  and  cooling  processes 
induced  by  picosecond  pulses  provide  a  unique  oppor¬ 
tunity  to  study  the  thermodynamics  of  transient  phase 
transitions  and  the  kinetics  of  crystal  regrowth  and 
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phase  transformation  under  extreme  conditions 

[1.2], 

Previously,  we  have  reported  that  at  higher  fluences. 
evaporation  of  material  from  the  irradiated  surface  can 
take  place.  The  morphologies  of  evaporation  or  re¬ 
moval  of  material  induced  by  a  20-ps  laser  pulse  under 
various  conditions  have  been  studied  [6].  When  a 
silicon  surface  covered  by  a  layer  of  Si02  is  irradiated 
with  a  picosecond  pulse  at  a  fluence  above  certain 
threshold  value,  the  oxide  layer  can  be  removed  after 
the  pulse.  The  sharp,  reproducible  threshold  laser 
fluences  to  induce  melting  of  the  silicon  surface, 
cracking  of  the  Si02  layer,  and  complete  removal  of 
SiOj,  for  samples,  with  oxide  layers  of  different 
thicknesses  are  reported  in  this  paper.  We  also  describe 
picosecond  time-resolved  reflectivity  measurements  on 
silicon  surfaces  covered  by  oxide  layers  of  known 
thickness  in  comparison  with  data  taken  on  bare 
silicon  surfaces  to  demonstrate  the  very  fast  processes 
of  heat  diffusion  and  material  evaporation  in  different 
laser  fluence  regimes. 


Experimental 

Silicon  wafers  of  (100)  and  ( 1 1 1 )  surfaces  were  used  in 
this  experiment.  However,  the  data  do  not  depend  on 
the  surface  orientation  except  that  the  amorphous 
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Fig.  1.  Reflectivity  of  a  bare  silicon  surface 
vs.  pump  fluence  at  various  probe  time 
delays.  Both  pump  and  probe  are  20-ps 
pulses  at  532-nm  wavelength 


patterns  induced  by  the  20-ps  pulses  at  /  =  532  nm  are 
easier  to  observe  on  a  ( 1 1 1 )  surface.  The  oxide  layers  of 
well-controlled  thickness  characterized  with  an  argon 
laser  [7]  were  grown  by  heat  treatment  in  dry  02,  at 
960r  or  U00°C,  respectively,  for  appropriate  time 
durations.  In  this  experiment,  we  exclusively  used  the 
20-ps  pulses  at  532-nm  wavelength  from  the  second 
harmonic  of  a  mode-locked  Nd :  YAG  laser.  The  time- 
resolved  reflectivity  measurements  were  performed 
with  the  same  experimental  setup  described  in  [4], 
Briefly,  the  sample  is  heated  by  a  pump  pulse  focused 
to  a  diameter  of  300  pm  at  the  (1/e)  intensity  contour  of 
its  Gaussian  spatial  profile.  The  probe  pulse  at  the 
same  wavelength  (/.  =  532  nm)  focused  to  a  diameter  of 
30  pm  at  the  exact  center  of  the  pump  spot  is  incident  at 
26”  to  the  normal  of  the  sample  surface.  The  variable 
time  delay  of  the  probe  pulse  with  respect  to  the  pump 
pulse  can  be  extended  from  -50  ps  to  2.5  ns. 

Results 

Here,  we  discuss  the  results  on  three  different  samples: 
bare  (111)  silicon,  (111)  silicon  with  a  thin  Si02  layer  of 
410 A  thickness,  and  (111)  silicon  with  a  thick  Si02 
layer  of2100  A.  Figure  1  shows  the  reflectivity  of  a  bare 
■silicon  surface  as  a  function  of  pump  fluence  at  various 
probe  delays  longer  than  100  ps  when  the  pump  and 
probe  pulses  are  temporally  completely  separated.  A 
sharp  rise  of  the  reflectivity  to  69%,  characteristic  of 
the  metallic  reflectivity  of  molten  silicon  at  an  incident 
angle  of  26°,  indicates  the  occurrence  of  melting  at  a 
threshold  fluence  of  0.2  J/cm2  [4, 5],  At  lower  fluences. 


the  reflectivity  increases  from  an  initial  value  of  31.5% 
gradually  with  increasing  pump  fluence  because  of 
heating  of  the  silicon  surface  below  the  melting  thres¬ 
hold.  At  higher  fluences,  the  decrease  of  the  reflectivity 
with  increasing  pump  fluence  at  100  and  200  ps  probe 
delays  indicates  overheating  of  the  molten  surface  with 
excess  fluence  above  the  melting  threshold.  This  re¬ 
duced  reflectivity  recovers  to  the  69%  level  as  the 
overheated  liquid  silicon  surface  cools  down.  This 
process  takes  about  300  ps  at  a  fluence  below  0.5  J  /cm 2 . 
Cooling  and  resolidification  of  the  molten  silicon 
surface  are  driven  by  diffusion  of  heat  into  the  bulk 
substrate.  At  fluences  just  above  the  melting  threshold, 
these  processes  may  take  place  right  after  the  pulse. 
When  the  resolidification  front  reaches  the  optical 
penetration  depth  of  molten  silicon  (about  90  A  at 
532  nm),  the  reflectivity  starts  to  drop.  This  explains 
the  gradual  drops  of  the  reflectivity  after  1  ns  probe 
delay  at  fluences  between  0.2  and  0.26  J/cm2.  Without 
the  precise  knowledge  of  the  initial  melt  depth  at 
various  pump  fluences,  we  cannot  derive  an  accurate 
value  of  the  upper  limit  of  the  regrowth  speed. 
However,  if  we  estimate  a  conservative  melt  depth  of 
400  A,  the  regrowth  speed  in  the  initial  stage  of 
resolidification  may  well  be  above  30  m/s.  This  high 
regrowth  speed  results  in  amorphization  of  the  silicon 
surface  observable  at  fluences  between  0.2  and 
0.26  J/cm2  [6, 8].  At  fluences  above  0.26  J/cm2  cooling 
takes  longer  times  as  indicated  by  the  long-lasting  high 
reflectivity.  The  slower  regrowth  speeds  then  allow  the 
liquid  surface  to  resolidify  into  crystalline  phase  at  a 
threshold  fluence  of  0.26  J/cm2. 
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Fig.  2.  Reflectivity  of  a  silicon  surface  with  a  SiO- 
layer  of4IO  A  thickness  vs.  pump  fluence  at  various 
probe  time  delays.  Both  pump  and  probe  arc  20-ps 
pulses  at  532-nm  wavelength 


When  a  silicon  surface  has  a  SiO,  layer,  the  optical 
reflectivity  can  be  changed  drastically  due  to  multiple 
interferences  in  the  oxide  layer  [7]  although  SiO,  does 
not  absorb  light  at  532  nm.  In  our  experiment,  the 
sample  with  a  thin  SiO,  layer  has  an  initial  reflectivity 
of  26.5%  at  532  nm  and  the  one  with  a  thick  Si02  layer 
has  an  initial  R  =  30%.  which  are  both  very  close  to  the 
bare  silicon  reflectivity  of  3 1 .5% .  For  the  sample  with  a 
thin  oxide  layer,  melting  also  occurs  at  a  threshold 
fluence  of  0.2  J  cm2  as  indicated  by  the  sharp  rise  in  the 
plot  of  reflectivity  vs  pump  fluence  at  various  probe 
delays  in  Fig.  2.  At  fluences  below  0.2  J/cm2  heating  of 
the  silicon  surface  also  increases  the  reflectivity  gradu¬ 
ally  with  increasing  pump  fluence.  However,  no  change 
in  the  behavior  of  multiple  interferences  as  a  function 
of  probe  time  delay  was  observable.  At  a  very  reprodu¬ 
cible  threshold  fluence  of  0.26  J/cm2,  we  started  to 
observe  cracking  and  partial  removal  of  the  Si02  layer. 
Complete,  clean  removal  of  the  Si02  layer  occurs  at  a 
threshold  of  0.3  J/cm2.  Because  Si02  softens  at  1853 
and  melts  at  1983  K  [9]  which  are  both  higher  than  the 
melting  point  of  silicon  at  1685  K  and  because  Si02 
does  not  absorb  light  at  532  nm,  at  fluences  between  0.2 
and  0.26  J/cm2  the  silicon  surface  melts  under  the  solid 
oxide  layer  and  resolidifies  into  amorphous  phase 
afterwards.  In  this  fluence  regime,  the  high  reflectivity 
level  rises  during  the  pulse  and  reaches  a  maximum 
value  of  61%  at  100  ps  delay,  which  is  lower  than  the 
reflectivity  of  the  bare  molten  silicon  surface  because  of 
the  multiple  interferences  in  the  oxide  layer.  However, 
the  reflectivity  does  not  stay  at  the  maximum  value,  but 
it  keeps  on  changing  and  starts  to  drop  after  lOOps 
when  heat  diffusion  from  the  molten  silicon  surface 


into  the  Si02  overlay  changes  the  optical  property  of 
SiO,  gradually.  This  drop  of  reflectivity  then  stops  at 
500  ps  and  the  reflectivity  stays  constant  with  time 
until  the  resolidification  front  under  the  molten  silicon 
layer  comes  within  one  optical  absorption  depth  from 
the  surface  after  about  1  ns,  as  discussed  previously. 
The  temperature  dependence  of  the  heat  difTusivity  in 
Si02  is  not  available  to  us.  If  we  take  the  only  available 
value  £)  =  0.006  cm2/s  at  room  temperature  [10].  we 
estimate  a  time  duration  on  the  order  of  1  ns  for  the 
heat  to  diffuse  across  the  410A-thick  SiO,  layer. 
Presumably,  the  time-resolved  reflectivity  data  show 
that  heat  diffusion  establishes  a  uniform  temperature 
profile  in  the  Si02  layer  within  500  ps.  Then  heat 
diffusion  into  the  bulk  silicon  substrate  is  solely 
responsible  for  further  cooling  of  the  hot  silicon  surface 
and  the  Si02  layer  afterwards. 

At  higher  fluences,  heating  of  the  Si02  layer  through 
heat  diffusion  from  the  overheated  molten  silicon 
surface  may  raise  the  temperature  to  the  softening 
point  or  even  the  melting  point  of  Si02.  Meanwhile, 
the  vapor  pressure  of  the  overheated  silicon  surface 
also  increases  with  increasing  pump  fluence.  The 
softened  or  molten  Si02  may  react  with  the  molten 
silicon  to  form  SiO  which  then  evaporates  [8].  At  even 
higher  fluences,  direct  evaporation  of  Si02  or  blow-off 
of  Si02  by  silicon  vapor  pressure  can  take  place.  After 
the  Si02  layer  is  completely  removed,  the  surface  then 
shows  the  high  reflectivity  R  =  69%  of  the  bare  molten 
silicon  surface.  The  reflectivity  plots  in  Fig.  2  then 
demonstrate  the  time  scale  of  the  evaporation 
processes.  At  fluences  between  0.26  and  0.3  J/cm2 
where  the  removal  of  Si02  is  never  complete  and 
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Fig.  3a-c.  Reflectivity  changes  as  a 
function  of  the  probe  time  delay  at  three 
different  sample  surfaces  irradiated 
with  a  20-ps.  532-nm  pulse  at  a  flucncc 
of 0.23  J /cm2: (a) bare  silicon  surface. (hi 
silicon  surface  with  a  SiO,  layer  of 
410  A  thickness,  (c)  silicon  surface  with 
a  SiO  2  layer  of  2 1 00  A  thickness 


dean,  the  reflectivity  recovers  gradually  to  higher 
levels  that  never  reach  69%.  In  this  fluence  regime  all 
the  processes  of  heating,  softening,  cracking,  and 
partial  removal  of  Si02  can  happen  and  the  changes 
in  the  reflectivity  are  very  complicated.  Above 
0.3  J/cm2  where  complete  removal  of  the  Si02  layer  is 
observed,  the  reflectivity  eventually  reaches  the  69% 
level  after  certain  time  delay  which  depends  on  the 
pump  fluence.  At  fluences  just  above  0.3  J/cm2  the 
evaporation  process  may  take  2  ns.  At  fluences  above 
0.4  j/cm 2  it  takes  150  to  200  ps  to  complete  the 
evaporation. 

The  data  on  the  sample  with  a  thick  Si02  layer  of 
2100  A  thickness  also  demonstrate  very  interesting 
features  of  the  heat  diffusion  and  evaporation  pro¬ 
cesses.  However,  because  of  the  larger  thickness  of  the 
Si02  layer,  multiple  interferences  and  heat  diffusion 
together  make  the  reflectivity  strongly  dependent  upon 
probe  time  delay.  The  data  became  too  complicated  for 
a  plot  as  in  Fig.  1  or  2.  In  this  sample,  the  threshold 
fluence  for  melting  is  increased  to  0.22  J/cm2  presum¬ 
ably  because  more  heat  is  diffused  into  the  thicker  Si02 
layer  during  laser  heating  of  the  silicon  surface.  At 
fluences  below  this  threshold,  the  reflectivity  increases 
gradually  with  increasing  pump  fluence.  In  contrast  to 
the  trace  in  Fig.  2  for  the  sample  with  a  thin  Si02  layer, 
the  reflectivity  now  changes  with  time  delay  even 
below  the  melting  threshold  because  the  2100  A  Si02 
layer  thickness  is  closer  to  the  probe  wavelength  of 
5320  A  so  that  changes  in  the  optical  property  in  this 
layer  due  to  heating  induce  observable  effects  from 
multiple  interferences.  At  higher  fluences,  removal  of 
the  Si02  layer  is  also  observed  with  reproducible 


thresholds  of  0.33  J/cm2  for  cracking  and  partial 
removal  and  0.39  J/cm2  for  complete  removal 
Although  the  layer  is  much  thicker,  the  process  of 
evaporation  and  removal  also  takes  only  a  few 
hundred  picoseconds  to  about  2  ns,  depending  on  the 
pump  fluence.  At  a  time  delay  of  1-1.5  ns,  the  reflectiv¬ 
ity  data  become  very  scattered  due  to  scattering  of  the 
probe  light  by  the  evaporated  Si02  particles  in  front  of 
the  sample  surface.  The  data  become  clean  again  after  a 
longer  delay  when  these  particles  are  completely 
removed. 

At  a  fluence  between  0.22  and  0.33  J/cm2,  the  reflectiv¬ 
ity  data  of  the  sample  with  a  thick  oxide  layer  show 
early  cooling  of  the  molten  silicon  surface  due  to 
diffusion  of  heat  into  the  Si02  layer  before  the  heat 
diffusion  into  the  bulk  silicon  substrate  starts  to  cool 
down  the  molten  silicon  layer.  This  phenomenon  is 
particularly  obvious  at  fluences  just  above  the  melting 
threshold  of  0.22 J/cm2. 

For  the  purpose  of  comparison,  the  reflectivity  changes 
of  the  three  samples  pumped  at  the  same  fluence  of 
0.23  J/cm2  are  shown  in  Fig.  3  as  a  function  of  probe 
time  delay.  The  reflectivity  of  the  bare  silicon  surface 
rises  abruptly  to  the  reflectivity  of  molten  silicon 
during  the  pulse,  stays  until  resolidification  front 
driven  by  heat  diffusion  into  the  bulk  substrate  reaches 
an  optical  depth  from  the  molten  surface  after  1  ns,  and 
then  drops  gradually  as  the  whole  molten  layer 
resolidifies.  For  the  sample  with  a  thin  Si02  layer,  the 
reflectivity  also  rises  during  the  pulse  as  the  silicon 
surface  melts  under  the  solid  Si02  layer.  It  reaches  a 
maximum  of  61%  at  lOOps  time  delay  and  falls  as  the 
Si02  layer  is  heated  by  the  heat  diffusion.  The  reflecti  v- 
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and  evaporation  processes.  Transparent  layers  with 
lower  melting  and  boiling  points  than  the  melting 
point  of  the  substrate  may  also  be  useful,  as  would  a 
time-resolved  mass-spectropraphic  studs  of  the 
evaporating  species. 
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Note  added  in  proof:  I.  In  order  to  avoid  interference  of  pump 
and  probe,  the  pump  beam  is  s-polanzed,  whereas  the  probe  is 
p-  polarized. 

2.  The  error  limits  on  the  R  values  are  ±  0.5%  at  zero  or  very 
low  pump  fluences  and  ±  1%  at  higher  pump  fluences  close  to 
or  above  the  melting  threshold. 


ity  then  stays  constant  from  500  ps  to  1  ns  time  delay 
after  the  oxide  layer  acquires  a  uniform  temperature 
profile  and  before  resolidification  of  the  underlying 
molten  silicon  layer  starts  to  reduce  the  reflectivity. 
The  reflectivity  trace  of  the  sample  with  a  thick  SiO_, 
layer  reveals  a  completely  different  process  of  cooling. 
It  rises  somewhat  more  gradually  to  its  maximum 
value  of  63%  at  100  ps  time  delay,  indicating  that  more 
heat  is  being  taken  away  by  diffusion  into  the  thick 
Si02  overlay  during  the  pulse.  After  it  reaches  the 
maximum  value,  it  drops  quickly  and  reaches  a  very 
low  value  after  500  ps.  Because  of  the  thickness  of  the 
oxide  layer,  a  uniform  temperature  profile  can  never  be 
established  in  the  layer  through  heat  diffusion  in  a  few 
nanoseconds.  This  thick  layer  also  has  enough  capac¬ 
ity  to  absorb  the  heat  released  by  the  molten  silicon 
layer  during  cooling  and  resolidifying.  Therefore,  this 
reflectivity  trace  shows  that  the  molten  silicon  layer  is 
now  cooled  and  resolidified  mainly  by  heat  diffusion 
into  the  thick  Si02  layer  in  less  than  500  ps.  rather  than 
being  cooled  by  heat  diffusion  into  the  bulk  substrate 
which  should  happen  later.  The  final  reflectivity  at  a 
level  of  33%  simply  demonstrates  the  reflectivity  of  a 
warm  amorphous  silicon  layer  covered  by  a  warm 
Si02  layer  with  a  nonuniform  temperature  profile, 
after  the  molten  silicon  layer  has  resolidified  into  the 
amorphous  phase. 

This  example  shows  that  the  behavior  of  coated 
surfaces  can  give  additional  information  about  heating 
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Using  a  novel  three-pulse  technique,  essential  information  about  the  density,  optical  effective 
mass,  and  kinetics  of  laser-generated  plasmas  in  silicon  has  been  obtained. 

PACS  numbers:  72.30.  -f  q,  79.20.Ds,  71.35.  +  z 


Numerous  investigations  of  the  change  in  optical  prop¬ 
erties  of  silicon  induced  by  a  strong  pump  or  “heating”  pulse 
have  recently  been  carried  out. 1-4  By  varying  the  wavelength 
of  the  probing  pulse  and  its  time  delay  with  respect  to  the 
picosecond  pump  pulse,  the  contributions  due  to  a  variation 
in  lattice  temperature  and  those  due  to  changes  in  carrier 
density  have  been  separately  determined.9 

The  real  part  of  the  dielectric  function  at  a  probing  fre¬ 
quency  a> 

*  "iini  1  ~  *f2(7’)l  —  (1) 

increases  with  the  lattice  temperature  T(dn/dT >  0,  tc2  <  1 ) 
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and  decreases  with  the  number  of  electron-hole  pairs  N.  As 
long  as  there  is  no  reliable  information  about  the  modifica¬ 
tion  of  the  optical  reduced  mass  m*  =  1  +  mV  V ' 

high  densities  and  carrier  temperatures  T,  available,  the  re¬ 
flectivity  data  of  laser-generated  plasmas  in  silicon  allow 
only  the  determination  of  N  /m*. 

The  imaginary  part  of  the  dielectric  function 


f" 


=  2«i(7>|.(7')  + 


4tre2N  /  1 

o5  V  m*,  <r,  > 


m*»  (t, 


y) 

(2) 


increases  with  Wand  T;  however,  as  a  third  unknown  param¬ 
eter,  the  scattering  times  of  the  carriers  (r,,  rj  averaged 
over  their  energy  distribution  come  into  play.  Thus  the  ex¬ 
perimental  determination  of  e’  and  e"  does  not  provide  suffi¬ 
cient  information  to  solve  separately  for  N,,  N,,  and  m*. 

In  this  letter  results  obtained  with  a  three-pulse  tech- 
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FIG  I  Time-resolved  reflectivity  and  transmission  in  SOS  induced  by  the 
first  pump  pulse  at  532  nm  with  $0  mJ/cm!  and  probed  at  the  same  wave¬ 
lengths  at  different  time  delays  {full  line).  Changes  induced  by  the  second 
p  imp  beam  at  1064  nm  with  25  mJ/cm1  (dotted  line) 


FIG  2  Relative  changes  of  reflectivity  and  transmission  probed  after  400 
ps,  induced  by  a  second  pulse  at  1064  nm  with  a  fluence  F‘„  =  40  mJ/cmv 
The  dff  and  dT  values  are  measured  as  a  function  of  its  delay  time  dr-,. 


nique  are  reported.  The  first  pulse,  with  a  fluence  of  50 
mJ/cm2  in  the  visible  (/l  =  530 nm)of20-ps duration,  creates 
the  carrier  plasma.  A  second  pulse  at  A  =  1.06  /rm  of  30-ps 
duration  adds  energy  to  the  plasma  by  the  free-carrier  ab¬ 
sorption  processes  without  changing  its  density.  A  third 
pulse  (A  =  530  nm)  probes  the  resulting  changes  in  reflectiv¬ 
ity  and  transmission  of  a  0.5-/zm-thick  silicon-on-sapphire 
(SOS)  sample.  Details  of  the  experimental  geometry  have 
been  given  in  Ref.  5. 

The  time  delay  between  the  first  two  pulses  is  kept  fixed, 
Jr2i  =  100  ps  Then  the  reflectivity  and  transmission  as  a 
function  of  df,„  the  time  delay  of  the  probe  pulse,  is  shown 
in  Fig.  1  for  two  cases.  The  drawn  lines  are  valid  in  the  ab¬ 
sence  of  the  second  pulse.  This  behavior  has  been  explained 
previously3  in  terms  of  decaying  plasma  density  and  with  the 
asymptotic  values  of  R  and  7" indicating  the  change  in  lattice 
temperature.  The  dotted  lines  are  valid  when  the  IR  plasma 
heating  pulse  of  25  mJ/cm2  is  applied.  The  data  give  the 
variations  dR  and  d  T,  due  to  the  change  in  lattice  tempera¬ 
ture  There  is  no  evidence  for  an  increase  in  N  /m*  due  to  the 


second  pulse.  The  data  indicate  that  the  energy  relaxation 
time  tw  of  the  carriers  with  the  lattice  is  faster  than  impact 
ionization. 

In  Fig.  2  the  probe  pulse  is  fixed  at  At}l  —  400  ps.  The 
change  in  reflectivity  AR  is  measured  as  a  function  ofd/2l.  It 
therefore  measures  the  variation  in  the  lattice  temperature. 
Clearly,  the  additional  heating  by  the  second  pulse  of  40  m)/ 
cm2  is  maximum  for  At2,  =  10  ps  and  decays  as  the  plasma 
density  created  by  the  first  pulse  decays. 

These  data  are  consistent  with  a  model  in  which  the 
energy  relaxation  time  tw  between  carriers  and  lattice  is 
taken  to  be  short  compared  to  the  pulse  duration.  In  fact,  we 
have  put  r  =  0.  Instantaneous  thermal  enhancement  of  the 
indirect  absorption  at  0.53  fim  is  assumed.  Thus  the  genera¬ 
tion  rate  of  electron-hole  pairs  experiences  a  significant  in¬ 
crease6  during  the  picosecond  excitation  at  0.53  #  With 
these  assumptions,  the  highest  possible  plasma  densities  in 
the  absence  of  impact  ionization  are  calculated.  The  results 
for  a  fluence  level  of  F =  100  mJ/cm2  are  shown  in  Fig.  3, 
where  the  plasma  density  at  the  surface  N,  and  the  surface 


-TOO  100  200  300 


TIME  (psac) 


FIG.  3.  Calculated  temporal  profile 
of  the  plasma  density  N  ■»  N,  =  N„ 
and  lattice  temperature  at  the  sur¬ 
face  of  silicon  during  and  following  a 
20  ps,  0. 1  i/cm1  pump  pulse  at  0. 53 

fim 
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temperature  Ts  are  plotted  versus  time.  The  plasma  densit> 
reaches  the  maximum  with  ~6.7  x  10^ 20  cm  - 3  a  few  pico¬ 
seconds  after  the  intensity  maximum  of  the  Gaussian  tempo¬ 
ral  excitation  profile  centered  at  At  =  0.  The  plasma  density 
is  strongly  reduced  by  Auger  recombination  (y  =  4x  10-3I|. 
The  rise  and  peak  of  the  lattice  temperature  are  retarded 
with  respect  to  the  plasma  features. 

Optical  probing  of  the  plasma  generated  at  this  fluence 
level  yields  an  upper  limit  of  X  /m*  values  with  X  /m* 
=  4x  1048  g_  1  cm-5.3,4,7  The  maximum  possible  value  of 
the  optical  mass  at  this  fluence  level  is  m*<0.18  m0.  Thus, 
the  data  are  consistent  with  a  nearly  constant  reduced  opti¬ 
cal  mass  m*  =  (m,‘  1  +  m,'  ')“ 1  =0.15  m0,  expected  for  a 
plasma  in  thermal  equilibrium  with  the  lattice. 

The  solid  lines  in  Fig.  2  are  calculated  on  the  basis  of 
this  model  with  the  free-carrier  absorption  cross  section  at 
1.06  nm  taken  to  be  a  —  2.3 X  10” 17  (7Y300).  This  coeffi¬ 
cient  is  four  times  larger  than  that  reported  at  lower  carrier 
concentrations.8  The  carrier  momentum  relaxation  times  r, 
and  rh,  occurring  in  Eq.  (2),  are  drastically  reduced  at  high 
temperatures  and  high  plasma  densities;  alternatively,  the 
Drude  model  does  not  adequately  describe  the  absorption.  A 
more  complete  theoretical  discussion  will  be  presented  else¬ 
where. 

In  conclusion,  the  three-pulse  experiments  support  the 


picture  that  the  carrier  plasma  rapidly  cools  ofT  by  transfer¬ 
ring  energy  to  the  lattice.  Plasma  heating  by  free-carrier  ab¬ 
sorption  does  not  lead  to  significant  impact  ionization.  The 
“simple  heating"  model  equations  yield  the  correct  plasma 
density  and  temperature  on  a  picosecond  time  scale.  The  real 
part  of  the  dielectric  constant  is  consistent  with  a  low  value 
of  the  effective  mass,  and  the  imaginary  part  indicates  a  shor¬ 
tening  of  the  momentum  relaxation  times. 
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We  have  observed  plasmon  resonances  with  1.9-^m  and  2.8-/xm  probe  pulses  in  silicon  and 
germanium,  excited  by  25-ps  pulses  up  to  40  mJ/cm2  at  0.53  and  1.06  jim,  respectively.  Firm 
values  of  N/m*  are  derived. 

PACS  numbers:  72.30.  +  q,  71.45.Gm,  52.25.Ps,  52.50.Jm 


The  kinetics  of  laser-generated  electron-hole  plasmas  in 
semiconductors  and  their  possible  influence  on  phase  transi¬ 
tions  have  been  of  considerable  interest,  particularly  in  the 
area  of  laser  annealing  of  ion-implanted  materials.1  One  of 
the  most  direct  techniques  which  can  be  used  to  obtain  the 
carrier  density  and  its  temporal  and  spatial  evolution  in¬ 
volves  time-resolved  free-carrier  infrared  spectroscopy2 
with  the  ambipolar  plasmas  probed  in  reflection  and/or 
transmission.  This  technique  has  been  used3  4  in  nanosecond 
and  picosecond  time  scales  to  observe  transient  plasmon  re¬ 
sonances  at  5.4  and  10.6 /xm  in  Ge,  Si,  and  GaAs  for  plasma 
densities  of  lO'MO20  cm-3.  To  investigate  higher  densities 
in  Si,  several  researchers  have  used  high  picosecond  or  fem¬ 
tosecond  excitations  with  probe  pulses  at  wavelengths  <  1 .06 
fj. m.5  Although  a  free-carrier  induced  drop  in  the  reflectivity 
is  observed  in  all  cases,  the  enhanced  reflectivity  associated 
with  the  plasmon  resonance  is  not.  Due  to  a  lack  of  knowl¬ 
edge  of  the  detailed  dielectric  function  at  high  excitation  lev¬ 
els,  it  is  difficult  to  extract  quantitative  information  about 
the  plasma  parameters  from  such  experiments.  However, 
even  assuming  the  validity  of  Drude-Zener  formalism,  opti¬ 
cal  measurements  allow  only  the  evaluation  of  N/m*,  where 
N  is  the  plasma  density  and  m*  the  electron-hole  reduced 
conductivity  effective  mass.6  From  previous  experiments 
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employing  a  1.06-/xm  picosecond  probe  a  normalized  den¬ 
sity  of  N/m*  =  3.4  X 1048  g-1  cm-3  has  been  found  for  a 
0.53-/xm,  25-ps,  100-mJ/cm2  excitation  pulse.6  According  to 
the  Drude  model  the  plasmon  resonance  would  therefore  be 
expected  at  a  wavelength  of 

A  =  c/e(m*/Nl>)ul(wem  )l/J~2/xm,  (1) 

where  =  1 1.8  is  the  high-frequency  dielectric  constant. 

In  this  letter  we  report  the  direct  observation  of  the 
plasmon  resonances  by  using  picosecond  probe  pulses  at  1 .9 
and  2.8  fim.  A  firm  value  of  N/m*  is  derived. 

A  passively  mode-locked  Nd:yttrium  aluminum  garnet 
laser  was  used  to  produce  single  30-ps,  1 .06-//m  pulses  which 
had  Gaussian  spatial  and  temporal  profiles.  These  pulses 
were  respectively  focused  to  spot  sizes  of  800-/xm  diameter 
(e~ 1  intensity  points)  on  ( 1 1 1 )  surfaces  of  nearly  intrinsic  Ge 
and  Si  wafers.  The  back  surface  of  these  samples  was  optical¬ 
ly  diffuse  and  so  prevented  undesired  interference  effects. 
The  generated  plasmas  were  probed  in  reflection  {R )  and  in 
transmission  (7")  using  1.9  and  pulses  which  were 

derived  from  a  portion  of  the  initial  1.06-/xm  pulse  by  the 
stimulated  Raman  scattering  in  50-cm-long  cells  containing 
50  atmospheres  of  H2  and  CH4,  respectively.  The  first  Stokes 
component  from  the  H2  gas  and  the  second  Stokes  compo¬ 
nent  from  the  CH4  gas  were  selected  using  interference 
filters.  The  pulse  widths,  although  not  measured,  are  esti¬ 
mated7  to  be  less  than  30  ps.  The  probe  pulses,  whose  spatial 
profile  was  measured  to  be  close  to  diffraction  limited,  were 
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Bulk  Silicon  wafer 


Bulk  GERuknium  WAFER 


FIG.  1.  Reflectivity  and  transmission  of  bulk  Si  at  1.9  pm  as  a  function  of 
probe  time  delay  after  a  0.53-pm  pump  pulse  at  0.04  J/cm:. 

focused  to  spot  sizes  of  <200  pm  at  the  center  of  the  excita¬ 
tion  region;  their  fluence  was  kept  less  than  1  mJ/cm2  to 
avoid  undesired  heating  effects.  Three  PbS  detectors  were 
used  to  measure  reflected,  transmitted,  and  reference  probe 
pulse  energies,  while  a  calibrated  vacuum  photodiode  was 
used  to  monitor  pump  energies.  The  signals  from  all  four 
detectors  were  handled  by  a  computer-automated  data  ac¬ 
quisition  system.  Standard  picosecond  pump  and  probe 
techniques  were  used  to  determine  R  and  T  as  a  function  of 
time. 

Figures  1-3  display  typical  results  for  the  time-resolved 
reflectivity  and  transmission  signatures  in  silicon  and  ger¬ 
manium  for  different  pump  fluences.  The  curves  accompa¬ 
nying  the  data  are  guides  to  the  eyes.  The  pump  fluence  was 
kept  below  the  threshold  fluence  for  melting  (200  mJ/cm2  in 
Si  at  0.53  pm,  75  mJ/cm2  in  Ge  at  1.06  pm)  to  avoid  any 
contributions  from  metallic  or  molten  regions.  Complete 
theoretical  interpretation  of  these  results  will  be  given  else¬ 
where.*  Here  we  wish  to  illustrate  the  salient  features  for 


FIG.  2.  Reflectivity  and  tmumudOfi  of  bulk  Si  at  2.8  pm  at  a  function  of 
probe  time  delay  after  a  0. 53pm  pump  pubed  at  (a)  0.0 1 5 J/em1  and  (b)  0.04 

l/em’. 
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FIG.  3.  Reflectivity  and  transmission  of  bulk  Ge  at  1.9  jxm  and  2.8  pm  as  a 
function  of  probe  time  delay  after  a  1 .06pm  pulse  of  fluence  SO  mJ/cm1 

different  excitations  and  probe  conditions.  Figure  1  shows 
the  reflectivity  R  at  1.9  pm  as  a  function  of  time  for  Si  for  a 
40-mJ/cm2  pump  fluence.  At  40  mJ/cm2,  and  below,  the 
transient  R  shows  a  broad  single  minimum.  Within  the 
Drude-Zener  model,  the  true  minimum  is  reached  when  the 
real  part  of  the  dielectric  constant  is  unity,  which  occurs  at  a 
critical  density  Ne  =  (1  —  l/e_  ]NP,  where  Nf  is  given  by 
Eq.  (1).  The  value  of  the  reflectivity  at  the  minimum, 
R  =  0.12,  is  determined  by  the  plasma  damping  mechanism 
and  the  temporal  width  of  the  probe  pulse.  The  nonzero 
width  of  the  probe  pulse  unfortunately  partially  washes  out 
the  actual  transient  optical  response  of  the  plasmon  reso¬ 
nance.  However,  at  least  for  the  broad  minimum  in  A,  at  low 
fluences,  these  effects  are  small.  Assuming  that  the  convolu¬ 
tion  effects  do  not  change  the  value  of  the  minimum  reflec¬ 
tivity,  a  value  of  N/m*  =  3.6x  104*  g_l  cm-3  is  derived.  If 
one  assumes  for  the  effective  mass,  m*  =  0. 123/no,  (the  low 
density  value)  where  m0  is  the  free-electron  mass,  the  plasma 
density  would  be  4X 1020  cm-3  at  the  pump  fluence  of  40 
mJ/cm2.  Theoretical  considerations9  indicate  that  m*,  and 
consequently  N,  may  be  up  to  factor  2  higher  at  this  level  of 
irradiation. 

The  transmission  data,  normalized  to  the  bulk  trans¬ 
mission  of  an  unexcited  sample,  provides  complementary 
information  to  that  of  the  reflectivity,  since  it  is  sensitive  to 
the  total  number  of  carriers  in  the  bulk.  The  normalized 
transmission  T  is  given  by 

r(f)  =  (0.7)-'[l  - i? (r )]exp(  -  erjr”2V(Jc)rf*),  (2) 

where  0.7  is  the  surface  transmission  of  the  unexcited  sam¬ 
ple.  This  equation  is  justified  if  <r,  the  free-carrier  absorption 
cross  section,  is  independent  of  depth  x.  The  transmission 
drops  to  a  minimum  of  10%  at  40  mJ/cm2.  This  transmis¬ 
sion  minimum  is  located  at  a  slightly  positive  time  delay 
between  the  excitation  and  the  probe  pulses.  This  is  a  direct 
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consequence  of  the  temporal  integration  (nonzero  probe 
pulse  widthl.  The  free-carrier  induced  transmission  change 
recovered  more  slowly  than  the  corresponding  reflectivity 
trace,  due  to  the  fact  that  the  plasma  density  at  the  surface  is 
higher  than  the  average  density  in  the  bulk. 

Figure  2  shows  typical  experimental  results  obtained  at 
2.8/zm,  where  the  plasmon  resonance  is  easily  observed.  The 
15-mJ/cm2  data  are  similar  to  that  of  Fig.  1  and  show  a 
broad  single  minimum.  At  40  mJ/cm2,  the  transient  R  shows 
two  shallow  minima  on  both  sides  of  a  peak  near  the  zero 
time  delay.  During  the  leading  and  the  trailing  edge  of  the 
pump  pulse,  the  carrier  surface  density  passes  twice  through 
the  critical  density  Nc .  Deconvolution  of  the  data  with  re¬ 
spect  to  the  nonzero  probe  pulse  width  indicates  that  the 
peak  reflectivity  is  >  0.9.  The  slower  decay  of  R  at  2.8  ^m  is 
as  expected,  since  the  reflectivity  at  the  longer  wavelength  is 
more  sensitive  to  a  given  density  of  carriers.  The  clear  experi¬ 
mental  observation  of  the  plasmon  resonance  at  2.8  ^m  indi¬ 
cates  that  the  N /m*  value  passes  twice  through  the  value  of 

1.8  X 104®  g“ 1  cm-3  on  each  side  of  the  maximum  value  of 
3.6X104®  g-,cm“3  at  the  same  fluence  of  40  mJ/cm2. 
Apart  from  discrepancy  near  zero  delay,  where  large  fluctu¬ 
ations  were  noted,  the  transmission  data  at  40  mJ/cm2  are 
consistent  with  Eq.  (2)  with  a  =  5x  10“ 7  cm2.10 

The  data  taken  with  the  Ge  sample  are  interesting,  since 
they  reveal  apparent  saturation  of  the  carrier  generation  rate 
which  has  been  the  subject  of  speculations  over  the  past  few 
years."-12  Figure  3  shows/?  and  T  data  at  1.9and  2.8/zm  for 
Ge  at  50  mJ/cm2, 30%  below  the  threshold  for  bulk  melting. 
With  a  pump  absorption  depth  of0.7  fim,  the  carrier  genera¬ 
tion  rate  at  1  06pm  is  higher  than  that  of  100  mJ/cm2  at  0.53 
pm  in  Si.  Also  the  Auger  recombination  coefficient13  of 
2  X  10“ 31  cmVs  is  a  factor  of  2  smaller  than  in  Si,  as  is  the 
reduced  effective  mass.  However,  as  shown,  the  plasmon  res¬ 
onance  is  not  observed  at  1 .9  pm  and  is  barely  discernible  at 

2.8  pm.  Theoretical  considerations  for  germanium  indicate 
that  the  reduced  effective  mass  is  not  more  than  20%  greater 
than  its  low  density  value  m*  =  0.08m0.  A  detailed  fit  of  the 
data  shows  that  the  maximum  density  obtained  is 
JV=  2.5x  1020  cm"3,  while  Np  =  2.0X  1020  cm'3  for  2.8 
pm.  The  observations  at  2.8  pm  are  consistent  with  those  at 

1.9  pm  where  no  reflectivity  peak  is  observed,  as  this  would 
require  a  density  at  the  surface  of  4.5  X  1020  cm-3.  The  low 
densities  are  indicative  of  intravalence  band  absorption1214 
or  saturation  of  the  valence-conduction  band  transition  due 
to  the  Burstein-Moss  effect.  A  model  which  neglects  these 
processes  would  lead  to  a  larger  carrier  density  at  the  surface 
(>  1021  cm-3),  which  is  not  observed.  The  low  minimum 
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value  of  the  transmission,  shown  here  only  for  1.9  pm,  also 
indicates  a  stronger  plasmon  damping,  due  to  the  added  con¬ 
tribution  of  intravalence  band  transitions.  Lastly,  it  should 
be  noted  that  the  melting  threshold  of  Ge  at  1.06  pm  is  con¬ 
sistent  with  the  absorption  coefficient  and  the  known  melt¬ 
ing  temperature  of  1210  K. 

In  summary,  we  have  shown  experimentally  that  the 
plasmon  resonances  are  located  between  1.9  and  2.8  pm  for 
picosecond  excitation  pulses  of  40  mJ/cm2  on  silicon  and  of 
50  mJ/cm2  on  germanium.  In  silicon  the  maximum  plasma 
density  at  the  surface  is  4X  1020  cm'3  if  variation  of  the 
effective  mass  is  negligible,  although  both  may  be  up  to  a 
factor  of  2  higher.  In  germanium  the  maximum  carrier  den¬ 
sity  is  about  2-3  X 1020  cm-3  before  melting  occurs.  These 
data  show  that  fast  infrared  transient  spectroscopy  of  Si  and 
Ge  is  a  valuable  technique. 

This  research  was  supported  by  the  Office  of  Naval  Re¬ 
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Picosecond  Photoemission  Study  of  Laser-Indnced  Phase  Transitions  in  Silicon 
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Pulsed  laser- induced  phase  transitions  in  semiconductors  are  character¬ 
ized  by  an  nltrafast  energy  transfer  from  the  generated  electron-hole  plasma 
to  the  phonons.  Time-resolved  optical  measurements  provide  ample  evidence 
that  the  thermal  equilibrium  between  carriers  and  phonons  is  established 
within  a  few  picoseconds  [1] .  However,  optical  techniques  mainly  explore 
the  evolution  of  plasma  density  and  lattice  temperature.  Direct  information 
about  the  energy  content  of  the  electron-hole  plasma  and  the  energy  loss  due 
to  phonon  interaction  cannot  be  obtained.  Picosecond  photoemission  provides 
a  direct  approach  to  this  problem. 

According  to  the  well-known  generalized  Fowler-DuBridge  theory  [2,3], 
the  electron  emission  current  is  composed  of  photoelectric  and  thermionic 
contributions.  Under  high  excitation  levels  multiphotoelectric  emission  is 
readily  observable.  Its  power  law  dependence  is  characterized  by  the  number 
of  photons  necessary  to  bridge  the  energy  gap  between  the  ground  state  and 
the  vacuum  level  in  intrinsic  silicon.  A  variation  of  the  photoelectric 
yield  is  expected  as  soon  as  the  incident  laser  pulse  induces  structural 
changes  on  the  surface  of  the  semiconductor.  The  thermionic  contribution  is 
extremely  sensitive  to  the  energy  content  of  a  thermal ized  carrier  distribu¬ 
tion  characterized  by  a  temperature  Tc  =  T#  =  T^.  Picosecond  excitation 
provides  a  stringent  test  for  anomalous  heating  in  which  the  carriers  are  no 
longer  in  equilibrium  with  the  lattice  phonon.  Due  to  the  small  specific 
heat  of  the  electron-hole  plasma,  Tc  should  be  much  higher  than  the  lattice 
temperature  T,  as  soon  as  the  energy  is  not  transferred  to  the  phonons 
instantaneously. 

In  this  contribution  we  report  on  picosecond  photoemission  of  silicon, 
with  photon  energies  corresponding  to  the  second  harmonic  (E  -  2.33  eV)  and 
fourth  harmonic  (E  =  4.66  eV)  of  a  Nd:YAG  laser  pulse.  The  photoelectric 
response  is  studied  as  a  function  of  laser  fluence  up  to  where  Ftb  is 

the  threshold  value  for  phase  transition.  This  work  extends  the  previous 
analysis  of  Liu  et  al.  [4]  to  laser  fluences  four  orders  of  magnitude  below 
F^,  to  study  explicitly  the  power  laws  of  photoemission  and  to  monitor  the 
occurrence  of  anomalous  heating  even  at  low  carrier  densities.  In  this  con¬ 
text  we  used  exclusively  4.66  eV  pulses  to  create  carriers  with  large  excess 
energies.  Specific  attention  has  been  paid  to  the  emission  of  positive 
ions,  which  starts  in  the  vicinity  of  F..  and  exhibits  an  extremely  non¬ 
linear  fluence  dependence.  It  is  clear  that  this  emission  is  intimately 
correlated  to  the  structural  change  of  the  surface  taking  place  at  Fth‘ 

The  experiments  were  performed  in  a  vacuum  of  5  x  10  Torr.  Both 
Si(lll)  and  Si(100)  have  been  irradiated.  A  vacuum  diode  configuration  was 
used  to  measure  the  electric  charges.  A  1  mm  diameter  wire  placed  at  2.5-3 
mm  from  the  sample  surface  acted  as  a  charge  collector.  Both  sample  and 
wire  were  electrically  insulated  from  the  body  of  the  vacuum  chamber  and 


could  be  independently  biased  up  to  +  4  kV.  Two  configurations  for  tbe 
measurements  are  possible.  One  can  measure  either  the  charges  collected  on 
the  wire  or  the  ones  escaping  from  the  sample.  Both  methods  have  been  shown 
to  give  similar  results. 

The  fluence  dependence  of  the  photoelectric  emission  can  be  generally 
divided  into  three  distinct  regimes.  In  the  case  of  UV-irradiation  (E  .  = 
4.66  eV) ,  a  superposition  of  linear  jnd  quadratic  effects  can  be  observed 
below  F(4w)  =  F  .  (4w)/15  ~  2  mJ/cni  (regime  I).  The  amount  of  collected 
charges  is  completely  independent  of  the  collector  potential.  Clearly,  at 
this  fluence  level  we  operate  in  the  saturation  regime,  colecting  all  elec¬ 
trons  emitted  from  the  surface.  The  value  of  the  linear  photoelectric  yield 
in  electrons  per  absorbed  photons,  Y.  =  8  i  10~  »  is  in  agreement  with  pub¬ 
lished  data  for  silicon  exposed  to  air  [6],  The  experimental  value  for  the 
quadratic  effect,  3  x  10  cm2  s  is  consistent  with  results  obtained  by 
Bensoussan  et  al .  in  picosecond  experiments  at  lower  photon  energies  [7]. 

In  regime  I  the  nonlinear  emission  dominates  the  electron  emission  over  a 
wide  range  of  fluences  (10  -10  J/cnr).  The  data  are  highly  reproducible 

and  the  quadratic  quantum  yield  is  completely  insensitive  to  surface  clean¬ 
ing  treatments. 

Above  F(4u>)  =  F  (4w)/15  the  collected  charge  signal  starts  to  be 
dependent  on  the  collector  potential.  We  enter  the  space  charge  limited 
regime  (II)  of  the  vacuum  diode,  where  the  applied  field  is  screened  by  a 
cloud  of  electrons  in  front  of  the  surface  [5].  The  nonlinear  emission  is 
masked  by  this  space  charge  effect.  The  electron  emission.,  is  limited  to  the 
leading  edge  of  the  picosecond  pulse.  At  F(4w)  =  25  mJ/cmZ,  where  the  first 
structural  changes  at  the  surface  are  observed,  the  electron  emission  exhi¬ 
bits  an  extremely  nonlinear  increase  in  this  regime  (III).  The  space  charge 
field  is  lowered  by  the  simultaneous  emission  of  positive  ions.  Above 
2Fth(4a>)  equal  amounts  of  electrons  and  positive  ions  can  be  collected  with 
the  +  4  kV  applied  on  the  wire.  The  space  charge  fields  are  completely  neu¬ 
tralized  above  this  fluence  level.  We  attribute  the  sharp  incrase  of  pho¬ 
toelectrons  and  the  simultaneous  apearance  of  poaitive  ions  to  a  structural 
change  of  the  surface.  The  threshold  for  photoemission  is  abruptly  changed 
to  lower  values  and  positive  ions  can  be  extracted  from  the  highly  disor¬ 
dered  phase  by  the  space  charge  field  developed  in  front  of  the  surface. 

This  picture  is  consistent  with  melting  of  the  surface.  Similar  behavior 


has  been  observed  in  the  case  of  VIS-irradiation  (E 


ph 


2.33  eV). 


The  photoelectric  emission  at  2.33  eV  is  determined  by  a  strong  two- 
quantum  process  and  weak  higher  order  contributions.  In  contrast  to  tbe 
4.66  eV  irradiation,  the  quadratic  quantum  yield  depends  strongly  on  the 
surface  condition,  as  expected  from  a  near-threshold  excitation.  Because  of 
the  lower  photoelectric  current  density  below  the  fluence  threshold  for 
phase  transition,  space  charge  effects  are  less  pronounced  than  in  the  case 
of  DV-irradiation. 


The  lower  quadratic  quantum  yield  at  2.33  eV  is  preferable  for  the 
study  of  possible  thermionic  contributions.  Application  of  the  Richardson- 
Dushmann  equation  in  the  regime  I  gives  an  upper  limit  of  2200  K  for  the 
carrier  temperature  up  to  fluence  levels  of  100  mJ/cm  .  At  this  VIS  fluence 
the  lattice  temperature  is  raised  to  ~  800  K  and  the  carrier  density  N  is 
limited  to  6  x  10  ,  as  the  optical  swasurements  show  [8].  This  maximum 


splitting  between  the  carrier  temperature  T  and  the  lattice  temperature 
confirms  that  most  of  the  energy  deposited  into  the  electron-hole  plasma  is 
transferred  to  the  phonons. 

This  research  was  supported  by  the  U.  S.  Office  of  Naval  Research  under 
contract  N00014-83K-0030 . 
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Abstract:  The  emission  of  second  harmonic  radiation  in  reflection  from 
crystalline  GaAs  irradiated  with  20  ps,  530  nm  laser  pulses  has  been  measured 
for  incident  laser  fluences  far  exceeding  the  threshold  fluence  for 

permanent  reflectivity  changes.  The  results  are  consistent  with  the 
occurrence  of  surface  melting  during  the  laser  pulse.  Detailed  analysis  of  the 
second  harmonic  signals  reveals  an  upper  limit  of  2  ps  for  the  structural 
transition  associated  with  the  melting  of  the  surface. 


The  intensity  of  second  harmonic  generation  (SHG)  from  non- 
centrosymmetric  crystals  such  as  GaAs  should  show  a  drastic  decrease  when  a 
transition  to  a  centrosymmetric  phase  occurs.  This  has  been  observed  with 
nanosecond  resolution  by  S.A.  Akhmanov  et  al.^  when  the  laser  fluence 
exceeds  a  threshold  value  for  melting. 

In  this  letter  we  present  the  results  of  SHG  study  on  GaAs  with 
picosecond  excitation.  In  simple  single  shot  experiments,  the  SHG 

dependence  of  the  incident  laser  fluence  is  monitored  under  appropriate 
geometrical  conditions.  The  deviation  from  the  quadratic  power  law  is 

analyzed  by  model  calculations  based  on  the  assumption  of  a  transition  to  a 
centrosymmetric  phase  as  soon  as  a  critical  fluence  during  the  20  ps  laser 
pulse  is  accumulated  on  the  surface  of  GaAs. 

In  our  analysis  we  use  mainly  (110)  surfaces  excited  with  single 

picosecond  pulses  at  532  nm  with  the  direction  of  the  electric  field  parallel  to 

a  <111>  crystal  direction.  The  second  harmonic  from  a  passively  mode-locked 

Nd:YAG  laser  system  is  focussed  on  the  GaAs  surface  at  an  angle  of  45°  of 

incidence.  The  pulse  duration  is  t  =  12.5  ps  at  the  1/e  points  of  the 

Gaussian  temporal  profile.  The  second  harmonic  signal  at  266  nm  is  detected 

in  the  direction  of  specular  reflection  by  a  photomultiplier  tube.  The 

reflected  exciting  radiation  is  suppressed  by  a  combination  of  filters.  Second 

harmonic  signals  down  to  =200  photons  are  detected.  The  inherent 

fluctuations  in  the  pulse  duration  of  the  laser  pulse  are  monitored  by 

(2) 

standard  .  A-techniques  Accidental  double  pulses  are  rejected.  Both 
precautions  yield  a  significant  smoothing  of  the  SHG  data. 


Generally,  two  separate  SHG  regimes  are  observed.  Below  the  critical 
2 

fluence  of  30  mJ/cm  at  532  nrn,  where  melting  and  subsequent  amorphization 

of  the  surface  occur,  the  signal  exhibits  the  expected  quadratic  power  law 

2 

dependence  I  =  n  I  ,  where  I  is  the  intensity  of  the  second  harmonic  pulse 

S  O  Q  S 

and  I  is  the  intensity  of  the  exciting  green  pulse.  Figure  1  shows  a  log-log 

plot  of  the  absolute  SHG  energy  emitted  from  a  (110)  surface  as  a  function  of 

the  excitation  fluence  F.  The  data  in  this  figure  are  integrated  signals  over 

space  and  time.  The  excitation  pulse  is  focussed  to  a  spot  size  of  200  ym  in 

diameter.  The  SHG  energy  increases  with  slope  2  up  to  a  fluence  level  of 

=  30  mJ/cm  .  The  SHG  efficiency  is  not  sensibly  affected  by  the 

generation  of  a  dense  electron-hole  plasma  or  by  lattice  heating  up  to  the 

melting  point.  A  value  of  ii  =  2.5x10  cm  /W  is  derived  from  the 

experiment.  Above  F^  the  SHG  signal  grows  more  slowly.  The  non- 

centrosymmetric  GaAs  structure  is  obviously  changed  during  the  excitation 

2 

pulse.  The  threshold  fluence  value  of  30  mJ/cm  is  in  excellent  agreement 
with  calculated  fluence  levels  necessary  to  melt  GaAs  surfaces^.  The  solid 
curve  in  Figure  1  represents  calculated  SHG  values  where  a  cutoff  in  the  SHG 
efficiency  is  assumed  as  soon  as  the  critical  fluence  F^  is  assumulated  during 
the  pulse.  The  integration  over  the  Gaussian  space  and  time  profiles  smooths 
the  sharp  drop  of  SHG  emission  considerably.  At  fluences  far  above  the 
threshold  value,  F^,  for  melting  the  UV  signal  increases  again.  This 
puzzling  behavior  has  been  carefully  examined  by  looking  at  the  angular  and 
time  dependence  of  the  signal.  For  pump  fluences  exceeding  =3F^  an 
additional  UV  signal  is  detected,  whose  fully  isotropic  angular  distribution  is 
in  marked  constrast  to  the  specular  SHG  beam.  This  radiation,  which  is 
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0 


► 
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accompanied  at  high  fluences  by  a  readily  observed  blue  spark  at  the 
surface,  exhibits  a  decay  time  of  =  lys,  while  the  duration  of  the  SHG  signal 
is  limited  by  the  time  constant  (10  ns)  of  the  photomultiplier  tube.  To 

analyze  the  growth  of  this  additional  signal,  the  emission  from  a  (100)  GaAs 
surface  is  studied,  where  the  polarization  of  the  excitation  beam  is  kept 
parallel  to  a  <001>  axis.  Under  these  conditions,  the  SHG  signal  is  zero  and 
the  background  UV  emission  can  be  measured  as  a  function  of  the  exciting 
532  nm  pulse  fluence.  As  shown  in  Figure  1,  this  new  emission  does  not 

affect  the  SHG  signal  at  fluences  below  5F^.  However,  due  to  the  extreme 
nonlinear  behavior,  it  is  responsible  for  the  unexpected  increase  of  the  UV 
emission  at  higher  fluences.  A  further  study  of  this  new  emission  process, 
which  is  closely  related  to  the  evaporation  and  plasma  formation  in  front  of 
the  surface  is  beyond  the  scope  of  this  letter.  Here,  we  restrict  our  interest 
to  the  time  scale  on  which  the  surface  structural  changes  on  GaAs  occur. 
For  this  purpose,  this  isotropic  contribution  at  high  pump  fluences  is  simply 
subtracted  from  the  total  UV  signal  on  (110)  surfaces. 

To  enhance  the  deviation  from  th'fe  quadratic  power  law  above  F^  the 

spatial  resolution  of  the  experimental  setup  is  improved  by  inserting  a 

diaphragm  in  the  path  of  the  specular  SHG  beam.  In  Figure  2,  the  energy 

dependence  of  the  SHG  emitted  by  the  central  portion  (120  ym  in  diameter)  of 

the  excited  area  on  the  surface  is  shown.  Clearly,  the  saturation  of  the  time 

integrated  SHG  signal  in  the  vicinity  of  F^  is  more  pronounced.  These  data 

are  now  compared  with  model  calculations,  where  the  transition  from  an 

ordered  non-centrosymmetric  structure  to  a  disordered  liquid  phase  is 

simulated  by  a  SHG  efficiency  decaying  with  a  time  constant  x  .  The  SHG 

s 

intensity  is  given  by: 


Is(r,t)  =  Ti(r,t)  I  ( r, t) x 


where 


u(r,t)  = 


for  t<t  (r) 
'  m 


Ho  exp  (-t-tm(r)  )  for  t  >  tm  (r).. 


At  each  point  r  on  the  surface,  n  starts  to  decay  as  soon  as  the  initial 


fluence  is  reached 


'8(r' 


t)  dt  =  F. 


The  solid  curve  in  Figure  2  corresponds  to  a  stepwise  truncation  of  the  SHG 
emission  by  setting  x  =  0.  Around  the  threshold  value  for  melting,  when  the 
phase  transition  still  occurs  on  the  trailing  edge  of  the  excitation  pulse,  the 
data  are  well  described  by  the  curve.  However,  at  higher  fluences,  when 
F^  is  reached  during  the  leading  part  of  the  excitation  pulse,  the  SHG  data 
points  are  slightly  above  the  calculated  curve. 

In  order  to  further  elucidate  the  temporal  behavior  of  the  transition  to 
the  disordered  phase  the  data  of  Figure  2  are  represented  in  a  linear  scale  in 
Figure  3  and  compared  with  the  model  calculations  with  different  time 

constants  t  .  Despite  the  scattering  of  the  data  points  the  solid  curves  of 

Figure  3  show  clearly  that  the  SHG  emission  drops  with  a  time  constant  xs  <  2 

ps  as  soon  as  the  energy  necessary  for  melting  is  locally  delivered.  In 


Figure  3,  the  calculation  for  t  =  5  ps  is  also  shown,  demonstrating  the 
sensitivity  of  the  SHG  emission  process  towards  a  slightly  delayed  phase 
transition. 

In  conclusion,  the  SHG  results  obtained  with  a  single  beam  technique 
indicate  that  an  ultrafast  transition  to  a  centrosymmetric  phase  occurs  in 
crystalline  GaAs  under  picosecond  laser  irradiation.  The  SHG  emission  follows 
the  usual  quadratic  power  law  up  to  the  threshold  fluence  for  meliing.  Above 
this  level  the  data  are  consistent  with  an  ultrafast  transition  to  a  disordered 
phase  within  2  ps. 
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Figure  1  Energy  of  the  specular  emission  at  26G  nm  vs  incident  laser  fluence 
(X  =  532  nm)  on  GaAs.  The  laser  spot  diameter  is  200  ym.  Closed 
dots  refer  to  (110)  surface,  laser  electric  field  parallel  to  the  [111] 
crystal  axis'.  Open  dots  refer  to  (100)  crystal  surface,  laser 
electric  field  parallel  to  the  [001]  crystal  axis.  The  solid  curve 
represents  calculated  SHG  energy  with  -  0  as  discussed  in  the 
text. 

Figure  2  Same  as  Figure  1,  but  limiting  the  collection  of  the  SHG  signal  from 
a  disk  of  120  ym  in  diameter  centered  on  the  incident  laser  spot. 

Figure  3  Plot  of  the  data  of  Figure  2  on  linear  scales.  The  solid  curves 

refer  to  calculations  with  the  values  of  t  shown. 
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Nonlinear  Photoemission  from  Picosecond  Irradiated  Silicon 

R.fl.flalvezzi ,  H.fCurz  and  N.BIoembergen 
Division  of  Rpplied  Sciences,  Harvard  University, Cambridge 
Nassachusetts  02138,  USR 


Three  distinctly  different  regimes  of  photoelectric  emission  are 
observed  over  a  aide  fluence  range  of  UU- laser  pulses  irradiating 
single-crystal  silicon  samples.  The  role  of  the  electron-hole 
plasma  in  the  nonlinear  photoemission  is  demonstrated  by  temporal 
correlation  measurements,  The  diffusion  properties  of  hot  carriers 
are  analyzed  by  investigating  the  influence  of  energy  transport  by 
hot  carrier  diffusion  on  the  fluence  threshold  for  melting  with  U*J 
photons . 


Pacs  numbers:  79.60, -i 


Pulsed  I  asen- i nduced  modifications  of  sem i condutor  surfaces 
find  growing  technical  interest  for  the  development  of  novel 
material  processing  technologies.  In  recent  times  much  experimental 
effort  has  been  devoted  to  studying  laser-induced  morphology 
changes  by  optical  techniques  [1].  The  formation  and  relaxation  of 
electron-hole  pairs,  the  kinetics  of  the  ambipolar  solid  state 
plasma  and  the  energy  transfer  to  the  lattice  phonons  has  been 
investigated  by  pump  and  probe  experiments  with  ps  and  fs  pulses 
[1 ,2,3,1] . 

Charged  particle  emission  measurements  provide  an  alternative 
approach  to  the  study  of  the  fundamental  process  at  the  surface  of 
laser  irradiated  material.  Contrary  to  optical  measurements,  where 
bulk  processes  are  dominant  and  the  interpretation  is  complicated 
by  spatial  gradients  of  the  parameters  to  be  studied,  here  one 
probes  the  outermost  layer  of  the  surface.  The  strict  power  law  of 
single-  and  multiple-quantum  excitation,  described  by  the 
generalized  Fowler-DuBr idge  theory  [5],  can  be  modified  by  thermal 
broadening  of  the  energy  distributions  of  the  carriers  in  the 
ground  or  excited  states,  fit  very  high  carrier  temperatures  the 
high  energy  tail  of  the  energy  distribution  reaches  the  vacuum 
level.,  and  pure  thermionic  electron  emission,  as  described  by  the 
generalized  Richardson-Dushman  equation..  occurs.  Nonlinear 
photoemission  of  silicon  has  been  observed  under  nanosecond  laser 
excitation  [6],  The  underlying  mechanisms,  however,  could  not  be 
clarified,  in  recent  picosecond  photoemission  experiments  the 
Richardson-Dushman  equation  has  been  used  to  determine  the  upper 
limit  of  the  electron  temperature  of  picosecond-irradiated  silicon 
at  laser  fluences  much  below  [7]  and  above  [8]  the  threshold  values 


3 


for  melting. 

In  this  paper  tue  present  for  the  first  time  a  plausible 
explanation  of  the  nonlinearity  of  photoemission  in  picosecond 
irradiated  silicon,  By  using  correlation  photoelectric  techniques, 
®e  haoe  been  able  to  reveal  the  dominant  role  of  the  electron-hole 
plasma  in  a  tuio-step,  single-photon  sequence.  In  addition,  by 
taking  advantage  of  the  extremely  small  penetration  depth  of  the 
laser  frequency  used  (a'  =  5  x  10“^  cm),  the  distance  L,  mhich  hot 
carriers  are  able  to  diffuse  before  they  transfer  their  energy  to 
the  lattice  phonons,  is  estimated  at  the  critical  laser  fluence  for 
melting. 

In  our  experiment  the  fourth  harmonic  (hp  *  4.66  eU)  of  a  30 
ps  Nd:VflG  laser  pulse  is  focused  on  a  crystalline  silicon  target 
placed  in  10“^  torr  vacuum.  The  emitted  charges  are  collected  by  a 
diode  configuration  as  described  in  [7,8],  ft  bias  voltage  of  1  kU 
is  applied  to  ensure  the  collection  of  the  total  amount  of  charges 
emitted.  The  dependance  of  the  collected  charge  density 
(Coulombs/cm41)  on  the  incident  laser  fluence  i3  shorn  in  Fig,  1. 

Three  distinctly  different  photoelectric  regimes  are  observed. 
In  regime  I  be  lorn  1  mJ/cm^  UU- laser  fluence,  a  superposition  of 
linear  and  quadratic  response  is  observed  independent  of  the 
applied  collector  voltage,  The  electron  current  density  J(r,t)  can 
be  expressed  in  this  regime  by 


J(r,t)  =  l(r,t)  e/hv  +  k2  (l(r,t)  e/hp)2 
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where  e  is  the  electron  charge,  hi>  the  laser  photon  energy,  l(r,t) 
the  space-  and  time-dependent  laser  intensity  on  the  sample 
surface.  Integration  of  (1)  ouer  a  space  and  time  gauss i an  profile 
results  in 

g/np2  =  k,  (e/hv)  F  +  2~3/2«‘1/2  k2  (e/hv)2  F2  (2) 

where  Q  is  the  total  emitted  charge,  p  is  the  1/e  radius  of  the 
laser  beam  on  the  surface,  r  its  1/e  time  duration  and  F  the  laser 
fluence  (J/cm2).  The  observed  linear  quantum  yield  of  8  t  2  x  10'^ 
electrons  per  absorbed  photon  is  in  fair  agreement  with  published 
data  of  silicon  exposed  to  air  [9],  The  value  of  the  coefficient 
for  the  quadratic  process,  k2  *  3.6  x  10'13  cit^s/Coul  is  consistent 

with  extrapolated  data  obtained  in  nanosecond  experiments  [6].  The 
dashed  line  of  Fig.  1  shows  the  calculated  charge  density  using 
these  yield  data.  In  regime  II  the  photoemission  deviates  from  the 
calculated  behavior  and  is  progressively  reduced  by  space  charge 
effects.  The  collected  charge  density  is  dependent  upon  the 
applied  voltage.  By  the  rapid  accumulation  of  electrons  in  front  of 
the  surface,  the  applied  field  is  screened  and  the  collector 
current  drops.  In  regime  III  the  electron  emission  exhibits  a 
strong  increase  again  as  soon  the  critical  laser  fluence  for 
melting,  Ft^  =  25  mJ/cm^  is  reached.  The  sharp  rise  in  the 

electron  emission  is  accompanied  by  the  emission  of  positive  ions, 
which  are  collected  when  the  bias  voltage  is  reversed.  The  space 
charge  field  which  inhibits  the  electron  emission  in  regime  II,  is 
neutralized  by  the  injection  of  positively  charged  particles.  Thus, 


5 


the  collection  of  electrons  above  the  threshold  for  melting  is 
governed  by  the  availability  of  positive  ions  to  remove  the  space 
charge  field.  Therefore,  nearly  the  same  amount  of  electrons  and 
ions  is  collected  at  higher  fluences  independently  of  the  magnitude 
of  the  app lied  field. 

In  order  to  explore  the  temporal  evolution  of  the 
photoemission  process,  we  used  a  time  correlation  technique  in 
which  two  consecutive  ps  pulses  are  focused  on  the  same  area  of  the 
sample.  The  synergetic  effects  on  the  photoemission  are  measured  as 

a  function  of  time  delay  x  between  the  two  pulses. 

To  begin  with,  this  double  pulse  technique  is  applied  for  a 
general  survey  of  regime  I  and  II.  In  Fig.  2  the  ratio  of  total 
collected  charge  Q  to  the  sum  of  laser  energy  E^  =  Ej  +  E2  is 

plotted  versus  the  sum  of  laser  fluences  Ft  =  F|  +  F2  (J/cm^).  The 

x  =  ••  condition  is  simulated  by  spatially  separating  the  two 
laser  spots.  The  solid  line  in  Fig.  2  illustrates  the  behavior 
under  this  condition.  Below  F  =  3  mJ/cm  the  quantum  efficiency 
increases  linearly,  as  expected  from  a  quadratic  process,  fit  7 
mJ/cffl^  the  maximum  efficiency  is  reached  and  at  higher  fluences 
(regime  II)  the  amount  of  collected  charges  drops  again. 

fis  soon  as  the  temporal  separation  of  the  two  laser  pulses  is 
less  then  100  ps,  noticeable  changes  in  the  fluence  response 
occur.  The  slope  of  the  linear  increase  below  3  »J/c*r  increases 
significantly.  The  collected  charge  maximum  moves  towards  smaller 
fluences,  and  the  drop  in  the  space  charge  regime  is  more 
pronounced.  The  largest  reduction  of  collection  efficiency  occurs 
in  this  regime  if  the  two  pulses  strike  the  surface  simultaneously 


( -i:  =  0 ) ,  In  this  case  the  total  charge  emitted  by  the  tuio  pulses 
is  less  than  the  sum  of  the  charges  of  tuio  separated  pulses, 
because  of  increased  suppression  by  the  space  charge  cloud  in  front 
of  the  surface.  If  the  tuio  pulses  hit  the  same  area,  but  are 
separated  by  a  time  interual  sufficiently  long  to  remoue  the  space 
charge  generated  by  the  first  pulse,  then  the  combined  charge 
emission  should  be  equal  to  that  of  tuio  spatially  separated  pul3e3. 
It  i3  estimated  from  our  sample  and  anode  field  geometry,  that  the 
electron  cloud  is  remoued  from  the  illuminated  spot  in  about  ID'10 
s.  The  data  in  Fig.  2  appear  to  confirm  these  space  charge  effects. 
Longer  lasting  space  charges  resulting  from  thermionic  emission  are 
negl igible. 

The  main  part  of  our  interest,  houieuer,  is  focused  on  the 
increase  of  quantum  efficiency  at  fluences  at  which  space  charge 
effects  can  be  neglected.  To  explore  the  origin  of  the  quadratic 
photoemission  process  in  regime  I,  ate  performed  a  complete  time- 
re3olued  study  of  the  correlation  of  the  collected  charges.  In  Fig. 
3  the  linear  slopes  of  Fig.  2  are  plotted  uersus  the  time  delay 

between  the  two  UU  pulses.  The  data  are  normalized  to  the  x  m  « 
condition.  Despite  the  large  experimental  error,  the  temporal 
correlation  behavior  clearly  indicates  that  the  first  pulse 
activates  the  surface  for  the  second.  This  optical  activation 
persists  for  times  much  longer  than  the  pulse  duration  and  is 
closely  correlated  to  the  dynamics  of  the  elctron-hole  plasma 
formed,  fit  a  fluence  level  of  I  mJ/cnr  no  significant  lattice 
heating  occurs,  and  the  maximum  carrier  density  generated  at  the 
surface  is  too  low  for  Auger  recombination  processes  to  become 


? 


dominant  in  the  density  equation.  The  penetration  depth  of  1.66  eU 
photons  (a1  =  5  x  10  cm)  is  much  smaller  than  the  diffusion 

j  n 

length  LQ  =  (Darp) '  ^  on  a  picosecond  time  scale.  Under  these 

conditions,  an  analytical  solution  of  the  density  at  the  surface 
N(0,t)  can  be  derived  using  solutions  of  the  surface  heating 
problem  [10].  Under  irradiation  mith  a  laser  pulse  with  gaussian 
temporal  profile,  the  density  t i me-dependance  is  given  by: 

c-b) 

N ( 0 , t )  -  -  x  W  r\(t/t  )  (3) 

P  P 

where  I m ( 1 -R)  is  the  reflection  corrected  incidence  irradiance, 
rp  the  l/e  half  width  of  the  loser  pulse,  DQ  the  ambipolar 
diffusion  coefficient  and  t\(t^p)  a  strictly  time-dependent  form 
factor,  which  is  developed  in  [9].  Uith  DQ  *  20  cm^/s,  a  maximum 

density  of  5  x  101S  cm"3  is  reached  after  0.55  x  under  the 
excitation  levels  used  in  this  correlation  experiment.  The  temporal 
decay  is  solely  determined  by  t|(t/tp)  .  The  striking  resemblance 
between  the  temporal  behavior  of  the  photoelectric  correlation 
signal  and  temporal  shape  of  n  leads  us  to  the  conclusion  that  the 
nonlinearity  of  the  photoemission  regime  I  is  caused  by  the 
presence  of  the  electron-hole  pairs.  It  is  basically  a  two-step 
single-photon  excitation  process,  where  in  the  first  step  electrons 
are  excited  via  direct  transitions  into  high-lying  conduction 
bands.  Rapid  intercarrier  collisions  ensure  thermal izat ion  of  the 
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carrier  distr i but  ion,  and  phonon  assisted  inter-  and  intravalley 
relaxation  processes  lead  to  the  accumulation  of  a  large  carrier 
density  in  the  vicinity  of  conduction  band  minima.  Because  of  the 
lack  of  sufficient  electron-hole  recombination  channels,  their 
accumulation  is  mainly  determined  by  the  relaxation  rate  from  the 
higher-lying  states  and  diffusion  into  lorn-density  regions.  From 
the  intermediate  excited  states,  as  well  as  from  the  final 
reservoir  in  the  conduction  band  minimum,  the  electrons  are 
reexcited  by  a  second  single-photon  absorption  process  to  final 
states  above  the  vacuum  level. 

fl  second  point  of  interest  in  this  paper  concerns  the 
diffusion  properties  of  hot  carriers  during  UU  excitation.  The  rise 
in  lattice  temperature  during  UU  excitation  at  higher  fluence 
levels  depends  sensitively  upon  the  distance  L  over  mhich  carriers 
have  diffused  into  the  bulk  before  substantial  energy  is 
transferred  to  the  phonons.  In  a  simple  approximation  this  distance 
L  can  be  treated  as  a  diffusion  length  [II] 


L  =  (  0Q  r  ),/2 


(la) 


mith  the  diffusion  coefficient  in  the  classical  model  given  by 


Da  * 


2  k  T  v  x. 

c  e  h 


*  * 

m  r,  ♦  m 

•  h  n  e 


(1b) 


In  these  equations  all  quantities  refer  to  hot  carriers,  In 
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(la)  ^  is  the  energy  relaxation  time  for  such  carriers,  and  in  eq. 

(lb)  x and  ^h  are  the  momentum  relaxation  times;  m  e  and  m  ^  are 
the  effective  carrier  masses,  and  T  is  the  temperature  of  the  hot 

C 

carriers,  fis  a  consequence  of  hot  carrier  di ffusion, the  heating 
rate  at  the  surface  by  UU  pulses  is  reduced  by  a  factor  (L  a  4 
I)-1,  mhere  a  is  the  UU  absorption  coefficient. 

In  numerical  model  calculations  the  influence  of  the  critical 
hot  carrier  diffusion  length  L  on  the  heating  rate  and  the  laser 
fluence  necessary  to  melt  the  surface  have  been  analyzed.  The  laser 
fluence  required  to  raise  the  surface  lattice  temperature  to  the 
melting  point  has  been  computed  as  a  function  of  L.  Comparison  with 
the  experimental  value  Fth  =  25  mJ/cm2  to  melt  silicon  with  UU 
photons  indicates  a  hot  carrier  diffusion  length  of  L  «  1  x  10-6  cm 
or  a  tia%  product  of  10" 12  cm2.  From  optical  investigations  it  is 

known  that  the  energy  relaxation  time  x  is  in  the  order  of  1  ps. 
fis  an  immediate  consequence  the  threshold  value  analysis  reveals  a 
rather  low  ambipolar  diffusion  coefficient  DQ  *  1  cm2/®  for  the 

hot  carriers  (comparable  to  that  in  metallic  conductors).  This 
should  be  compared  with  the  room  temperature  value  at  low  densities 
(  0Q  =  20  cn^/s).  In  agreement  with  transport  measurements  at  high 

carrier  densities  [12],  but  in  disagreement  with  a  series  of 
theoretical  predictions  [11],  the  ambipolar  diffusion  of  hot 
carriers  is  lower  than  the  room  temperature  value,  ft  drastic 


decrease  of  the  momentum  relaxation  times  x  down  to  1  x  ICf15  s 

e,h 

occurs  in  UU  laser- induced  electron-hole  plasmas  at  fluences  close 
to  the  melting  threshold,  outweighing  the  linear  dependance  of  Dq 

on  the  electron  temperature  T  .  This  momentum  relaxation  time  is 

two  orders  of  magnitude  lower  then  the  ualue  determined  by  carrier 
-lattice  scattering  (  *  2  x  10-13  s).  Similar  low  ualues  have  been 
found  in  optical  invest i gat  ions  of  the  free  carrier  absorption 
cross  sections  at  high  excitation  level  [13].  These  extremely  short 
momentum  relaxation  times  are  consistent  with  the  fast  energy 
losses  of  the  electron-hole  plasma  to  the  phonons  observed  in 
earlier  investigations  using  2.33  eU  excitation.  Uhether  this 
reduction  of  the  hot  carrier  diffusion  mobility  is  due  to  the 
extremely  high  excess  energy  of  the  initially-created  electron-hole 
pairs  or  to  a  density-dependent  phenomenon  like  electron  hole 
scattering  cannot  be  answered  conclusively  at  the  present  state  of 
the  investigation,  flore  theoretical  work  on  the  basis  of  the  data 
presented  in  this  paper  is  necessary  to  clarify  this  point. 

In  conclusion,  three  distinct  regimes  in  the  photoemission  are 
observed  in  ps  experiments  using  4.66  eU  photons.  The  quadratic 
photoelectric  response  is  closely  related  to  the  formation  of  an 
electron-hole  plasma  during  the  pulse  providing  the  intermediate 
states  from  where  electrons  are  more  easily  excited  into  the  vacuum 
level.  The  distance  which  the  initially  hot  carriers  can  travel 
before  they  lose  their  energy  to  the  lattice  phonons  indicates 
extremely  short  momentum  relaxation  times  under  excitation  levels 


1 1 

close  to  the  threshold  fluence  for  melting. 
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Figure  Captions 


Fig,  1  Emitted  charge  density  versus  fluence  at  266  nm  for  S  i  ( 1 1 1 ) 
samples.  The  dashed  line  refers  to  the  calculated  linear 
and  quadratic  photoelectron  effect  with  the  coefficients 
given  in  the  text , 

Fig.  2  Ratio  Q/En  of  the  negative  charge  emitted  by  S i ( 1 1 1 )  over 
total  laser  energy  at  266  nm  versus  total  laser  fluence 
for  two  pulses  separated  by  the  interval  x  shown. 

Fig.  3  Correlated  photoelectron  signal  versus  daly  time  x  for  266 
nm  illumination  of  Si  (1 11)  in  region  I.  The  data  are 
normalized  to  the  Xs"  results. 
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PICOSECOND  PULSED  LASER-INDUCED  MELTING  AND  GLASS  FORMATION  IN  METALS 


Chien-Jung  LIN,  Frars  SPAEPEN  and  David  TURNBULL 

Division  of  Applied  Sciences,  Harvard  University,  Cambridge,  Mass.  02138* 


Irradiation  of  a  solid  surface  by  a  picosecond  laser  pulse  leads  to  melting  of 
a  thin  overlay,  and  subsequent  quenching  at  rates  of  10^‘K/s  or  higher.  This 
technique  is  also  one  of  the  few  available  for  superheating  a  crystal  above  its 
melting  temperature.  The  ultra-fast  quenching  rates  permit  the  formation  of 
new  metallic  glasses.  Work  on  the  Fe-B  and  Ni-Nb  system  is  reviewed,  and  new 
results  on  the  Mo-Ni,  Mo-Co  and  Nb-Co  systems  are  reported.  In  all  cases, 
glasses  are  formed  below  the  T0-line,  indicating  that  partitionless  crystalli¬ 
zation  is  usually  at  least  partially  diffusion  controlled  due  to  the  change  in 
short  range  order  upon  crystallization;  only  in  the  simplest  structures,  such 
as  disordered  f.c.c.  solutions,  are  the  partitionless  crystallization  kinetics 
fast  enough  to  preclude  glass  formation. 


1.  INTRODUCTION 

When  a  solid  metallic  surface  is  irradiated  by  a  laser  pulse  of  duration  tD, 

1  y 
the  energy  is  deposited  in  a  layer  of  thickness  a  ,  where  q  is  the  absorption 

coefficient;  for  metals,  a"'  is  on  the  order  of  1 00& .  During  the  pulse,  a 
layer  of  thickness  (2Dthtp)  is  heated  ud  by  thermal  conduction,  where  Dtp  is 
the  thermal  diffusivity.  For  metals  Dph  is  on  the  order  of  0.5cm?/sec,  which, 
for  a  30ps  pulse,  corresponds  to  an  initial  heated  layer  of  about  500A  thick. 

For  a  typical  laser  fluence  of  about  0.5  J/'cm  and  a  metallic  reflectivity,  the 
temperature  rise  in  this  (molten)  layer  is  estimated  to  be  several  thousand 
degrees.  The  fraction  of  the  energy  lost  due  to  evaporation  is  small  ( <  1 0. ) 
since  it  is  limited  by  the  evaporation  kinetics  and  the  short  duration  of  the 
process.  The  rest  of  the  energy  is  conducted  into  the  bulk,  leading  to  further 

O 

melting  to  a  total  depth  on  the  order  of  1000A.  Further  cooling  then  leads  to 
resolidification  of  this  layer,  either  by  crystallization  or  glass  formation. 
Heat  flow  calculations  and  direct  observations  of  melting  (e.g.  in  Si,  using 
time-resolved  reflectivity)  show  that  the  total  lifetime  of  the  melt  for  a  30ps 
pulse  is  about  Ins;'  given  that  the  temperature  drop  is  on  the  order  of  lO^K, 

this  corresponds  to  an  average  cooling  rate  during  solidification  on  the  order 
1 2 

of  10  K/s.  Since  most  "conventional"  quenching  techniques,  such  as  melt  spin¬ 
ning,  produce  cooling  rates  of  about  10®K/s,  it  is  expected  that  pulsed  laser 
quenching  can  lead  to  the  formation  of  new  metastable  phases.  A  review  of  the 
fundamentals  of  pulsed  laser  interactions  with  solid  surfaces  has  been  given  by 

•This  work  has  been  supported  by  the  Office  of  Naval  Research  under  Contract  No. 
N00014-83-K-030. 
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?.  MELTING  AND  SOLIDIFICATION  KINETICS 

The  general  formalism  of  these  transformations  has  been  reviewed  by  Spaenen 
and  Turnbull.  '  The  velocity  of  the  crystal-melt  interface,  u,  may  be  approx¬ 
imated  by: 

u“u0(VV/Ti  «’! 

where  Tm  is  the  melting  temperature  and  T.  the  temperature  at  the  interface;  u0 
is  a  velocity  determined  by  the  kinetics  of  the  melting  or  crystallization  pro¬ 
cess:  for  pure  metals,  when  these  processes  are  collision-limited,  u0  is 
approximately  eaual  to  the  speed  of  sound,  u;;  for  alloys,  where  diffusive 
rearrangements  mav  be  important  u  «  D/1.  (D:  diffusivity;  1:  interatomic  dis¬ 
tance)  . 

The  heat  flow  associated  with  melting  or  solidification  leads  to  a  second 
expression  for  u: 

u  =  •<  V  trTMFc  (2) 

with  <:  thermal  conductivity;  V:  molar  volume;  f.T:  temperature  gradient  at 
tne  interface;  ,GH  :  molar  heat  of  crystallization.  The  combination  of  eqs .  (1) 
and  (?)  then  leads  to  a  determination  of  the  interface  temperature,  T^ .  If  the 
interface  kinetics  are  fast  compared  to  the  rate  of  heat  removal  (or  heat  sup¬ 
ply,  for  melting),  as  is  the  case  in  most  conventional  metal  processing,  the 
process  is  called  heat  flcr :  lirrirei,  and  T^  is  only  slightly  different  from  T^. 
If,  however,  the  interface  kinetics  are  slow  (e.g.  if  limited  by  diffusional 
rearrangements),  and  the  rate  of  heat  removal  (or  supply)  is  fast  (e.g.  due  to 
the  very  steep  thermal  gradients  in  ps  pulsed  laser  annealing),  the  process 
becomes  iiterface  limited,  and  Tj  is  very  different  from  T^. 

If  the  kinetic  processes  at  the  interface  are  collision-limited,  as  in  cry- 
sta  1  1  izat  ion  of  pure  metals,  T.  differs  from,  T^  by  only  a  few  tens  of  degrees 
even  in  ps  pulsed  laser  quenching.  It  should  be  noted  that  the  thermal  gra¬ 
dients  during  the  melt-in  phase  following  the  pulse  are  much  steeper  than  dur¬ 
ing  solidification:  V^T  lo\/(2Dthtp)1,/2  •»  lO^K/m.  According  to  eq .  (2) 
this  corresponds  to  a  heat-flow  limited  velocity  on  the  order  of  2000m/s,  i.e., 
on  the  same  order  of  u$ .  Eq.  (1)  therefore  predicts  a  substantial  overheat  at 
the  interface  Consequently,  picosecond  pulsed  laser  heating  is  one  of  the  few 
methods  available  for  superheating  metallic  crystals. 

lor  metallic  glass  formation,  it  is  clear  that  T^  must  fall  far  below  TM 
during  cooling,  and  that  hence  diffusion-limited  interface  kinetics  seem 
required.  The  solidification  velocity,  u,  in  ps  laser  quenching  can  be 
estimated,  from  the  melt  depth  ( -.  1000^)  and  the  lifetime  of  the  melt  ( ~1  ns), 


(  J  l  in  Ci  al  T’lt  t  ul,J  fKil\t  J  lj'<  r-lllj;a  »  J  nit 
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to  tp  at  lea'.t  10>/«.  The  distance  tnol  a'  alio/  etc,"  car  travel  diffj'ivel, 
during  the  passing  of  the  solidification  front,  If-'ui  ’£,  is  then  less  than  an 
interatomic  distance.  The  solidification  process  is  therefore  a  partition- 
less"  one.  Partitionless  crystall  1  tat ior,  is  only  tnermodyna-ica  1  ly  possible  if 
the  temperature  drops  below  the  so-called  "T^-line"  in  Fig.  1,  i.e.,  the  line 
where  the  free  energy  of  the  crystal  and  the  undercooled  liquid  solutions  are 
equa  1 .  J 

Compos-  »ion 


Cf  yslollizofton 

FIGURE  1 

Schematic  illustration  of  the  T0-lines  'equal  free  energy 
for  liquid  and  crystal)  in  a  simple  eutectic  phase  diagram. 

This  "solute  trapping"  process  has  been  discussed  in  detail  by  Aziz.6  If  the 
configurational  freezing  temperature  ("Tg")  of  an  alloy  lies  above  the  TQ 
lines,  glass  formation  is  expected  if  no  intermetal 1 ic  phases  can  be  formed; 
the  complex  structure  of  most  of  these  phases,  however,  usually  requires 
diffusion-like  rearrangements  that  slow  down  the  crystallization  kinetics. 

Since  partitionless  crystal  1 ization  into  the  primary  solutions  is  often  a  very 
fast  process,  it  has  been  proposed  that  the  T  -lines  represent  the  limits  of 

7  0  0 

glass  formation."'  However,  as  will  be  demonstrated  below,  this  is  only  the 
case  for  the  simplest  crystal  structures,  such  as  disordered  f.c.c.  alloys, 
where  the  kinetics  may  approach  the  collision  limit.  If,  as  is  often  the  case, 
reconstruction  of  the  short-range  order  is  required,  the  crystall ization  kinet¬ 
ics  are  at  least  partially  diffusion-limited  and  glass  formation  can  occur  even 
below  the  T  -1  ines  . 

3.  EXPERIMENTAL  PROCEDURE 

The  main  experimental  problem  in  studying  the  solidification  of  binary 
alloys  in  the  ps  pulse  regime  is  obtaining  a  homogeneous  melt  in  its  1  ns  life- 
time,  since  the  mixing  length  corresponding  to  this  time,  (Dt)  '  ,  is  only  30A. 
The  phases  in  the  alloy  must  therefore  be  dispersed  on  this  scale.  We  have 
solved  this  problem  ’  by  preparing  the  starting  alloy  as  a  ~1000A  thick  com- 
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positionally  modulated  film  with  a  wavelength  shorter  than  the  mixing  length. 
After  irradiation,  the  film  is  floated  off  the  substrate  and  is  ready  for 
transmission  electron  microscopy  without  further  thinning.  This  procedure  is 
also  very  flexible  in  that  it  allows  preparation  of  a  continuous  range  of 
compositions  using  only  two  sources  or  sputtering  targets. 


4.  RESULTS 

In  earlier  papers^''1^  we  showed  that  Fe-B  alloys  with  a  B  content  as  low  as 
6  at.'  can  be  made  glassy  by  psec  laser  quenching,  whereas  "conventional" 
quenching  requires  a  minimum  of  12  at.tB.  Below  5  at.tB,  the  alloys  solidify 
as  supersaturated  b.c.c.  crystals.  Estimates  of  the  TQ-line  in  the  Fe-B  sys¬ 
tem,  based  on  regular  solution  modeling  and  on  the  devitrification  mechanism, 
show  that  all  the  ne-j  alloys  obtained  by  laser  quenching  were  formed  below  the 
r  -1 ine. 

o  l  1 

Similar  observations  were  made  for  the  Ni-Nb  system,  where  glasses  were 

formed  by  laser  quenching  in  the  23-82  at.SNi  range,  with  (supersaturated) 

f.c.c.  and  b.c.c.  crystalline  solid  solutions  outside  this  range;  the  glass 
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FIGURE  2 

Molybdenum-nickel  phase  digram,  indi¬ 
cating  estimated  T0-lines  for  the 
primary  solid  solutions.  Glasses  have 
been  formed  by  laser  irradiation  at  the 
three  compositions  indicated  by  arrows. 


1  ? 

ing  is  only  40-70  at.SNi,  and 
corresponds  roughly  to  the  range 
between  the  TQ-lines.  Again,  laser 
quenching  demonstrated  that  glasses 
could  still  be  formed  far  beta:  the 
TQ-lines. 

We  also  report  here  new  results  on 
glasses  formed  in  the  Ni-Mo, Co-Mo  and 
Co-Nb  systems.  Figure  2  shows  the 
Mo-Ni  phase  diagrams  with  estimated 
To-lines  for  the  primary  solid  solu¬ 
tions.  For  three  compositions,  with 
30,  50  and  60  at.*Ni,  ps  laser 
quenching  produced  the  amorphous 
phase.  Note  that  in  two  cases  glasses 
are  formed  below  the  T0-line,  and  that 
the  <5 (MoNi )  intermetal  lie  is  not 
formed  upon  quenching  the  alloy  be¬ 
tween  the  T0-lines.  Similar  alloys, 
with  35,  50  and  65  at.fcNi,  had  been 
made  amorphous  by  liu  et  al}^  by  ion 
mixing  of  multi-layer  films  of  the 
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(a) 


FIGURE  4 


(b) 


Electron  diffraction  patterns  of  an  as-deposited  (a)  and 
laser-irradiated  (b)  compositionally  modulated  film  of 
average  composition  Nb^COgg. 


crystalline  elements.  In  the  60  at.tNi  alloy,  we  observed  some  microcrystals , 
but  only  near  the  center  of  the  laser-irradiated  spot;  in  the  30  at.tNi  alloy 
some  crystals  (b.c.c.)  were  observed  at  the  edge  of  the  spot,  but  not  in  the 
center;  in  the  50  at . ‘Ni  alloy  the  entire  spot  was  amorphous.  A  possible 
explanation  of  the  difference  in  the  occurrence  of  the  crystals  may  be  that 
crystals  at  the  edge  in  the  30  at.%Ni  alloys  grow  epitaxially  on  the  only 
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partly  melted  modulated  film  (the  beam  has  a  Gaussian  intensity  profile), 
whereas  in  the  60  at.Mii  alloy  the  higher  quench  rate  at  the  edge  suppresses 
growth  while  at  the  slower  cooling  center  new  crystals  can  nucleate  on  the 
substrate.  Figures  3  and  4  show  the  diffraction  patterns  of  as-deposited  and 
irradiated  Mo^Co^  and  Nb^pCo^Q  alloys.  In  both  cases,  fully  glassy  phases 
are  formed  upon  irradiation.  To  our  knowledge,  Co-Mo  amorphous  alloys  had  so 
far  not  been  produced  by  liquid  quenching;  Liu  et  at.  have  produced  amorphous 
alloys  with  65  and  35  at.JCo  by  ion  mixing.  It  is  interesting  that  the  as- 
deposited  modulated  Nb-Co  film  is  entirely  amorphous,  whereas  the  as-deposited 
Mo-Co  film  is  entirely  crystalline,  although  both  were  prepared  by  sequential 
sputter  deposition  of  the  crystalline  elements.  Since  both  systems  are  clearly 
glass  formers,  and  since  the  sputter-induced  mixing  is  probably  similar  for 
both  films,  the  occurrence  of  the  as-deposited  amorphous  film  may  be  the  result 
of  enhanced  solid  state  interdiffusion  between  the  layers,  as  observed  in 
La-Au.14 
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Resume  -  La  vitrification,  par  trempe  au  laser  a  duree  de  picosecondes , 
des  alliages  Fe-B  contenant  au  moins  de  5  at . %  B  est  expliquee  quanti- 
tativement  par  une  transition  d'une  croissance  cristalline  limitee  par 
collision  atomique,  a  une  croissance  limitee  par  diffusion  atomique.  la 
ligne  TQ  pour  la  solidification  de  la  phase  c.c.  dans  le  liquide,  sans 
redistribution  des  constituents,  est  calculee  utilisant  la  theorie  des 
solutions  regulieres  et  le  diagramme  de  phases. 

Abstract  -  Glass  formation  in  picosecond  laser  quenching  of  Fe-B  alloys 
containing  a  minimum  of  5  at.%  B  is  explained  quantitatively  by  a  transi¬ 
tion  from  collision-limited  to  diffusion-limited  crystal  growth.  The  T 
line  for  partitionless  growth  of  the  fee  phase  from  the  liquid  is 
calculated  from  regular  solution  theory  and  the  equilibrium  phase  diagram. 


1  -  INTRODUCTION 

The  process  of  energy  deposition,  heating  and  colling  in  pulsed  laser  processing  has 
been  discussed  in  detail  by  Bloembergen  [1]  and  others  [2-4] .  Figure  1  illustrates 


I  Energy  deposition 


2  Melt-in  3.  Cryttol  regrowth  or  glots  tor  mot  ion 
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1  ■  30- 100  pwc 


«■  lOOps-Int 


Fig.  1  -  Schematic  diagrams  illustrating  the  mechanism  of  pulsed  laser  quenching. 

the  main  stages  of  this  process  in  a  typical  picosecond  pulsed  laser  experiment: 

(i)  energy  deposition  in  a  layer  on  the  order  of  the  absorption  depth  (200  X) ,  which 
spreads  by  heat  conduction  over  a  layer  of  thickness  (2D  t  )1//2  of  about  500  % 


(Dty,:  thermal  diffusivity;  tp :  laser  pulse  length);  since  energy  losses  due  to 
evaporation  are  negligible  14 J ,  this  layer  melts  and  its  temperature  rises  to  several 
thousand  degrees;  (ii)  during  the  melt-in  period,  this  overheat  is  spent  as  heat  of 
fusion  of  more  substrate  crystal  to  a  total  thickness  of  about  1000  &;  due  to  the 
very  high  thermal  gradients,  the  estimated  melt-in  velocity  is  extremely  large 
(~1000  m/s) ;  (iii)  when  the  maximum  melt  depth  is  reached,  the  crystal  melt  inter¬ 
face  reverses  direction;  depending  on  the  relative  rates  of  crystal  growth  and  heat 
removal,  the  molten  layer  solidifies  as  a  regrown  crystal  or  a  glass.  The  competi¬ 
tion  between  the  latter  processes  are  analyzed  in  detail  in  this  paper  for  the  Fe-B 
system. 

2  -  THERMAL  PARAMETERS 

Although  a  number  of  sophisticated  analyses  of  the  heat  flow  during  laser  quenching 
are  available,  the  order  of  magnitude  of  most  of  the  relevant  quantities  can  easily 
be  estimated  from  essentially  dimensional  arguments.  They  are  listed  in  Table  I, 
together  with  the  corresponding  values  for  melt  spinning.  The  starting  point  is 

Table  I 

THERMAL  PARAMETERS  IN  MELT  QUENCHING 


Laser  Quenching 

Melt  Spinning 

Melt  temperature 

T  (K) 
m 

io3 

-7 

io3 

-5 

Melt  thickness 

d  (m) 

10 

5  *  io 

Temperature  gradient 

VT  (K/m) 

1010 

2  x  107 

Cooling  rate 

T (K/s) 

io12 

4  xio6 

Melt  lifetime 

T  (S) 

io-9 

not  applicable 

Isotherm  velocity 

uT (m/s ) 

100 

0.2 

Heat-flow  limited  crystal 

growth  velocity 

(m/s ) 

230 

0.5 

that  for  both  processes  the  melt  temperature,  T  ,  is  on  the  order  of  10  K,  but  the 
melt  thickness,  d,  in  laser  quenching  is  much  less  (1000  A)  than  in  melt  spinning 
(50  pm) .  The  temperature  gradient  for  heat  removal  can  then  be  estimated  as 
VT=Tm/d,  and  the  cooling  rate  as  T  =  D^VT/d  (Dth:  thermal  diffusivity,  about 
-5  2 

10  m  /s  for  iron  its  alloys).  The  lifetime  of  the  melt  is  then  T  =Tm/T;  this 
quantity  is  only  relevant  for  laser  quenching,  and  has  in  fact  been  checked  experi¬ 


mentally  by  time  resolved  reflectivity  measurements  on  Si,  which  has  thermal 
characteristics  similar  to  those  of  the  metals  [5,6].  The  velocity  of  the  isotherms 
in  the  specimen  (the  quantity  to  be  compared  to  the  crystal  growth  velocity  later) 
is  uT  =T/VT. 

3  -  THERMODYNAMICS  AND  KINETICS  OF  CRYSTAL  GROWTH 


The  crystal  growth  velocity  is  determined  by  the  balance  of  the  rer  oval  of  the  heat 
of  crystallization,  Ah,  and  the  kinetics  of  the  atomic  processes  at  the  interface. 
The  heat  flow  limited  crystal  growth  velocity  is; 


cal  values  fcr 


(x :  thermal  conductivity; 
Fe ,  usinc  AH  =15  kJ/mole, 


V  molar  volume).  Table  I 

=  50  W,  m  •  K ,  and  V  =  7 . 1  *  1 0” ^ 


lists  t 
m^/mole . 


The  crystal  growth  velocity  depends  on  the  interface  kinetics  and  the  thermodynamic 
driving  force  as  follows  [2]  : 


u  = 


(2) 


where  AG  is  the  free  energy  difference  between  crystal  and  melt,  and  the  prefactor, 
uQ,  depends  on  the  nature  of  the  atomic  rearrangement  needed  for  incorporation  of  a 
liquid  atom  into  the  crystal.  If  only  simple  collision  of  liquid  atoms  with  the 
crystal  surface  is  required  (as  in  pure  metals) ,  the  atomic  jump  frequency  is  the 
vibrational  frequency  in  the  liquid,  so  that  uD  in  this  case  is  close  to  the  speed  of 
sound,  us,  in  the  liquid.  In  the  analysis  below,  a  rough,  estimated  value  of 
3500  m/s,  independent  of  temperature,  is  used  for  us.  If  incorporation  of  the  atoms 
in  the  crystal  requires  diffusional  jumps  (i.e.,  changes  of  nearest  neighbors),  uD 
is  approximately  equal  to  D/X  (D:  diffusivity  in  the  liquid;  X;  jump  distance,  taken 
as  2.6  *10_1°  m  in  Fe-based  alloys).  Note  that  diffusion- limited  crystal  growth 
does  not  necessarily  imply  long-range  atomic  transport;  changes  in  short-range  order 
that  require  changes  in  coordination  can  also  only  occur  by  diffusional  jumps.  The 
diffusivity  used  in  the  analysis  below  has  the  Fulcher-Vogel  form,  typical  of  liquid 
metal  alloys: 

D  =  Do  exp(‘  (3) 

with  D  =1.4  * 10~8  m2/s,  T'  = 581K,  and  B  =  1300K. 
o  o 

The  value  of  P  is  a  rough  estimate,  similar  to  that  observed  in  measurements  on  other 
metallic  glasses  and  consistent  with  the  free  volume  model  17] .  Efe,  and  T^  were  then 
calculated  by  assuming  D=5  *10-9  m2/s  at  the  melting  point  of  iron  (1807K) ,  and 
D=10”X6  m2/s  at  the  crystallization  temperature  of  B-rich  Fe-B  glasses  (650K)  [8]. 


Since  the  isotherm  speed  is  several  100  m/s  in  the  ps  laser  quench,  the  crystal 
growth  velocity,  u,  must  be  at  least  this  high  to  prevent  glass  formation.  The  time 
required  to  crystallize  a  monolayer  at  this  speed,  X/u,  is  less  than  1  ps.  The 
distance  an  atom  can  diffuse  in  that  amount  of  time,  (DX/u)1/2,  is  then  less  than  an 
interatomic  distance.  Crystallization  must  therefore  be  partitionless ,  i.e.,  without 
change  in  composition.  In  alloys,  partitionless  crystallization  at  a  particular  com¬ 
position  is  thermodynamically  possible  only  below  a  temperature,  Tc,  at  which  the 
free  energies  of  melt  and  crystal  are  equal. 

Assuming  a  simple  regular  solution  model  for  both  melt  and  crystal,  the  driving  free 
energy  for  partitionless  crystallization  of  an  Fe-B  alloy  melt  to  a  supersaturated 
b.c.c.  phase  is 

AG  =  Gc  =  xFeAb°e  +xBAp°  +*FexB(ec  -V  (4) 


(Gc,  G^ :  free  energies  of  crystal  and  liquid,  respectively;  xFg,  xg:  atom  fractions 
of  Fe  and  B;  regular  solution  interaction  parameters). 

The  standard  free  energy  difference  can  be  approximated  by: 


AS,  .  (T  -T  .  ) 
f,i  M,  l 


(5) 


(AS,  entropy  of  fusion,  and  T„  melting  temperature  of  element  i).  The  T 

1,1  n,i  u 

temperature  is  found  by  setting  AG  equal  to  zero,  which  gives: 


o 


I 

Fe  E 


'c'V^FeVre^f  ,Fe4XBTK,E^£f  ,B 


x  As,  _  +x  As, 

Fe  f ,  Fe  B  f ,  B 


(6) 


The  interaction  parameters  can  be  determined  by  equating  the  chemical  potentials  of 

Fe  and  E  at  the  equilibrium  liquid  and  crystal  compositions,  which  can  be  found  or. 

the  phase  diagram.  For  the  Fe-E  systems,  Cc  -  C£  was  calculated  to  be  -92  kj/mole, 

independent  of  temperature.  The  other  parameters  can  be  found  in  all  standard 

references  [9]:  AS  =7.6  J/K  mole;  AS  =21.8  J/K  mole;  T  -  1807K; 

I / Fg  i f  B  M / Fc 

T  = 2300K.  The  result  is  shown  in  Fig.  2.  Note  that  only  the  b.c.c.  phase  is 
K ,  B 


Fig.  2  -  The  Fe-B  phase  diagrar. 
The  Tc  line  corresponds  to  the 
liquid  -*t  .c.c.  transition.  The 
contours  represent  the  crystal 
growth  speed  (in  m/s)  as  a  func¬ 
tion  of  composition  and  inter¬ 
face  temperature. 


being  -onsidered  here;  the  f.c.c.  Fe  phase  does  not  nucleate  or  grow.  Inserting 
Eq .  1-3)  into  Eq.  (4)  gives  an  expression  for  the  driving  force. 


AG 


‘V'Sf.Fe  <T 


T  ) 
o 


(7) 


Some  authors  110.11J  nave  assumed  that  partitionless  crystallization  is  always  a 
very  fast  process,  and  that  therefore  glass  formation  below  the  T0  line  is  prohi¬ 
bited.  The  Fe-B  system  is  one  of  many  in  which  this  has  experimentally  been  dis¬ 
proved 


■4  SXrLRi MENTAL  RESULTS 

Fe-:-  glasses  have  been  prepared  by  melt  spinning  in  the  12-28  at.%B  range  [8,12,13]. 
Bel.:  -i  i.2  at.  a  b.c.c.  supersaturated  solution  is  formed.  It  is  interesting  to 
note  th.it  upon  heating  glasses  with  less  than  18  at.%B  first  crystallize  to  the 
ai«t?stabie  h  c.c.  solution,  which  then  converts  to  the  stable  phase  mixture  of  a-Fe 
and  [81.  Our  calculation  shows  that  the  composition  corresponding  to  To=660K 

(the  tystailiZation  temperature  for  an  18  at.%B  glass)  is  16  at.%B.  Considering 
the  .approximate  natur?  of  the  regular  solution  models,  this  agreement  is  quite 
sa*-  i  s factory 

Clu.  ,  fox iiiat .on  in  the  Fe-B  system  has  been  studied  under  ps  laser  quenching  condi¬ 
tion  - ,  using  the  following  experimental  technique  [14,15).  On  a  copper  substrate, 
firs':  coated  with  a  1  |i  thick  A1  film,  a  1000  &  thick  compositionally  modulated 
film  consisting  of  alternate  layers  of  Fe  and  Fe^,B.,  was  deposited  by  dual  target 

/o  *4 


o 

DC  sputtering.  The  repeat  length  of  the  modulation  was  about  20  A.  By  varying  tne 
relative  thicknesses  of  the  Fe  and  layer,  films  with  average  compositions 

between  0  and  24  at.%B  were  prepared.  These  samples  were  irradiated  with  a  30  ps 
Nd : Y AG  (X  =1.06  pm)  laser  pulse,  with  a  beam  diameter  of  about  100  pm  and  an  average 
fluence  of  about  0.8  J/cm^.  The  composition  modulation  wavelength  was  chosen  on  the 
order  of  the  mixing  length  (DT)l/2  in  the  liquid  during  the  lifetime  of  the  melt 
(T  =  1  ns.  Table  I),  to  ensure  homogeneity.  After  dissolution  of  the  A1  film  in 
dilute  NaOH,  the  irradiated  film  could  be  investigated  directly  by  transmission 
electron  microscopy  without  further  thinning. 

These  experiments  showed  that  ps  laser  quenching  results  in  glass  formation  in  alloys 
containing  at  least  5  at.%B.  Alloys  containing  less  than  5  at.%B  crystallized  as  a 
b.c.c.  solid  solution.  The  solidification  morphology  of  the  4  at.%B  contained  a 
small  amount  of  glass/crystalline  mixture  near  the  top  surface  of  the  film  116] ; 
this  could  have  been  the  result  of  an  interfacial  instability  that  became  possible 
at  the  last  stages  of  solidification  due  to  the  decrease  in  the  crystal  growth 
velocity  [17];  the  associated  redistribution  of  B  then  could  result  in  glass  forma¬ 
tion  in  the  B-rich  regions. 


5  -  ANALYSIS 

The  sharp  transition  from  partitionless  crystal  growth  to  glass  formation  at  5  at.%B 
can  most  easily  be  explained  by  a  transition  from  a  collision-controlled  to  a 
diffusion- controlled  crystal  growth  mechanism.  The  short-range  order  in  a  metal- 
metalloid  glass  such  as  Fe-B,  and  hence  also  in  the  very  undercooled  melt,  is  known 
to  be  quite  strong:  each  metalloid  is  surrounded  by  metal  atoms  only,  in  a  coordi¬ 
nation  shell  similar  to  that  found  in  the  intermetallic  compounds  [18,19] .  The 
short-range  order  around  the  B  at~m  in  the  b.c.c.  solution,  either  substitutional  or 
interstitial,  is  clearly  very  different.  Reconstruction  of  the  cluster  containing 
the  B  atom  therefore  requires  diffusive  jumps.  The  formation  of  a  glass  with  a 
minimum  of  5  at.%B  suggests  that  the  size  of  this  cluster  is  20  atoms  (i.e.,  each  B 
atom  and  19  of  its  nearest  neighbors,  which  corresponds  to  the  first  and  about  half 
the  second  coordination  shell  of  the  B-atom  [20]  ).  The  Fe  atoms  outside  these 
clusters  are  assumed  to  make  collisional  jumps.  The  kinetic  factor  in  Eq. (2)  for  the 

crystal  grcwth  speed  can  therefore  be  written  for  x  ^0.05: 

B 

u  =  (20  x  )  y  +  (1-20  x  )u  .  (8) 

O  Jd  A  B  S 

For  xB^0.05,  the  crystal  growth  speed  is  equal  to  the  diffusive  speed  D/X. 

Figure  2  shows  contours  of  the  crystal  growth  speed,  u,  as  a  function  of  interface 
temperature  and  composition.  The  abrupt  transition  at  5  at.%B  is  obvious.  Figure  3 
shows  as  a  function  of  composition,  the  total  distance  grown  by  the  crystal  if  cooled 
at  a  constant  rate  form  the  liquidus  temperature  to  600K.  The  scale  on  the  left 
hand  side  is  for  the  "master  curve"  representing  the  product  of  this  distance  and 
the  cooling  rate.  On  the  right  hand  side,  the  actual  growth  distances  in  melt 
spinning  and  laser  quenching  are  shown.  At  5  at.%B,  the  growth  distance  in  ps  laser 
quenching  is  on  the  order  of  10  X;  most  of  a  1000  X  thick  layer  is  therefore  trans¬ 
formed  to  a  glass,  as  observed  in  our  experiments.  At  12  at.%B,  the  growth  distance 
in  melt  spinning  is  on  the  order  of  1  ps,  consistent  with  the  experimentally 
observed  glass  formation  in  50  pm  thick  ribbons.  Figure  4  illustrates  the  same 
points  in  a  slightly  different  way.  For  pure  Fe  and  the  4  at.%B  alloy,  the  crystal 
growth  velocity  rises  above  the  isotherm  velocity  in  ps  laser  quenching,  thus  pre¬ 
venting  glass  formation.  For  5  at.%B,  the  crystal  growth  velocity  rises  above  the 
isotherm  velocity  in  melt  spinning,  but  remains  below  the  one  in  laser  quenching. 

This  is  consistent  with  glass  formation  at  this  composition  under  the  latter  condi¬ 
tions,  and  failure  to  form  a  glass  under  the  former.  The  crystal  growth  velocity 
for  the  12  at.%B  alloy  lies  below  the  isotherm  velocity  for  melt  spinning,  in  agree¬ 
ment  with  the  observations. 
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Fig.  3  "  Crystal  growth  distance 
upon  continuous  cooling  of  ar. 

Fe-B  alloy,  as  a  function  of  com¬ 
position  (see  text  for  explana¬ 
tion  of  scales). 
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Fig.  A  -  Crystal  grcuth  velocity  as  a  function  of  interface  temperature  for  four 
compositions  in  the  Fe-B  system. 


There  exists  experimental  support  for  the  very  large  growth  velocities  for  pure  Fe 
shown  on  Fig.  4.  In  ps  laser  annealing  experiments  on  pure  Fe ,  the  original 
sputtered  film  had  a  grain  size  of  500  A.  After  irradiation  the  grain  size  was  1  pr 
or  more  115,21].  Since  these  grains  had  to  grow  within  the  1  ns  lifetime  of  the 
melt,  the  lateral  growth  velocity  must  have  been  at  least  500  m/s;  the  velocity  may 
in  fact  be  considerably  larger  yet,  since  the  lateral  growth  is  limited  to  the  very 
early  stages  of  the  crystallization.  Coriell  and  Turnbull  [22]  also  found  a  very 
large  kinetic  factor,  on  the  order  of  the  speed  of  sound,  in  their  analysis  of 
dendritic  growth  velocity  measurements  in  very  undercooled  Ni  melts. 


6  -  CONCLUS I ONS 

The  ps  laser  quenching  experiments  on  Fe-B  alloys  have  demonstrated  that  metallic 
glasses  can  be  formed  far  below  the  Tc-line.  This  has  also  been  observed  in  the 
Ni-Nb,  Mo-Ni  and  Mo-Co  systems  [3].  Partitionless  crystallization  is  therefore  not 
necessarily  a  fast  mechanism  that  forestalls  glass  formation  as  some  authors  have 
claimed  [10,11].  For  the  Fe-B  system,  the  glass  formation  range  can  be  qualitatively 
accounted  for  by  a  transition  from  collision-controlled  to  diffusion-controlled 
crystallization  at  5  at.%B.  The  analysis  gives  sensible  values  for  the  glass  forma¬ 
tion  range  in  melt  spinning  as  well.  Detailed  quantitative  agreement  on  the  latter 
process  cannot  be  claimed,  however,  since  the  results  at  low  temperature  depend 
rather  sensitively  on  the  choice  of  the  parameters  in  the  Fulcher-Vogel  equation  for 
the  diffusion  coefficient.  Crystal  nucleation,  which  is  probably  moire  important  in 
melt  spinning  than  in  laser  quenching,  has  not  been  considered  here.  Unless  the 
conditions  are  right  for  copious  homogeneous  nucleation  (including  the  transient 
effects  [23]),  taking  into  account  nucleation  effects  will  result  in  a  decrease  of 
the  fraction  crystallized. 
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1.  INTRODUCTION 

Compoaitionally  modulated  materials  are  fabricated  by  repeated  alternate 
deposition  of  very  thin  layers  of  different  composition.  Materials  with 
composition  modulation  repeat  lengths  between  a  few  %  and  about  100A  are  of 
interest  for  a  number  of  scientific  and  technical  applications  as 
semiconductor  devices,  as  z-ray  mirrors,  as  samples  for  the  study  of  diffusion 
at  low  temperature  [1],  and  as  starting  materials  for  alloying  in  the  liquid 
state  following  very  fast  pulsed  laser  heating  [2].  Some  of  their  physical 
properties,  such  as  their  elastic  moduli  [3]  and  magnetic  properties  M],  are 
unusual . 

These  materials  have  been  produced  by  a  variety  of  methods:  molecular 
beam  epitaxy  (5]  (mainly  for  semiconductors),  vapor  deposition  [6,7],  D.C.  and 
R.F.  sputtering  [8,9],  and  chemical  vapor  deposition  [10].  In  this  paper,  we 
report  on  a  new  method  that  employs  ion  beam  sputtering.  Its  main  advantages, 
as  will  be  explained  below,  are  its  simplicity,  its  flexibility  (e.g.  metals, 
semiconductors  and  insulators  can  be  deposited  under  identical  conditions), 
and  convenient  control. 
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2.1  Overall  design 

The  apparatus  consists  of  a  vacuo*  chamber  containing  two  ion  guns,  one 
aimed  at  the  target,  the  other  aiaed  at  the  substrate.  The  target  assembly 
contains  fonr  different  targets  that  can  be  aelected  externally  by  a  stepping 
■otor.  The  deposition  rate  is  monitored  by  a  crystal  oscillator.  A  feedback 
aystea  froa  the  deposition  aonitor  to  the  target  selection  stepping  notor 
controls  the  sequential  deposition  for  the  fabrication  of  compositional ly 
aodulated  aaterials. 

2.2  Chamber 

The  chamber  is  a  vertical  cylinder,  45cm  tall  and  41cm  in  diameter, 
made  entirely  of  stainless  steel.  Figure  1  shows  a  view  from  the  top  of  the 
cylinder.  The  bottom  plate  is  connected  to  the  vacuum  system.  The  top  plate 
contains  the  mechanical  feedthroughs  for  the  shutters  and  the  target  assembly, 
and  the  electrical  feedthroughs  for  the  ion  probes  in  the  gun  shutters. 

The  substrate  is  placed  on  a  water  cooled  support  block  on  a  removable 
flange.  This  flange  also  supports  the  crystal  sensor  and  a  second  substrate 
holder;  the  latter  is  a  cylinder  that  allows  external  insertion  of  a  heating 
element.  The  centers  of  the  targets  of  the  substrate  and  of  the  substrate 
flange,  and  the  ion  beams  are  all  in  the  same  plane.  The  centers  of  the 
crystal  sensor  or  the  hot  substrate  can  be  placed  in  the  same  plane  by 
rotating  the  substrate  flange  by  72°  or  144°  (2  or  4  boltholes),  respectively. 
The  substrate  flange  also  contains  a  shutter  by  which  the  substrate  or  the 
crystal  can  be  exposed,  either  separately  or  together. 


The  target  assembly  consists  of  a  hollow  4-inch  stainless  steel  cube 


attached  to  a  2-inch  vertical  stainless  steel  tube  that  goes  through  a 
Ferrof lnidics  SC2000C  hollow  rotary  feedthrough  in  the  top  plate  of  the 
chaaber.  The  targets  are  4-inch  square  plates  that  are  screwed  to  the 
vertical  faces  of  the  cube.  Cooling  water  is  brought  in  through  a  1-inch 
atainless  steel  tube  inside  the  2-inch  tube  and  is  removed  through  the  space 
between  the  tubes.  Outside  the  chaaber,  flexible  hoses  are  used  to  conduct 
the  cooling  water,  thus  allowing  360°C  rotation  of  the  target  asseably.  The 
target  asseably  is  surrounded  by  a  shield  that  ainiaizes  the  spread  of 
sputtered  aaterial  around  the  chaaber. 

All  the  ahutters  in  the  chaaber  are  aoved  through  Ferrof luidics 
SB250AN086  1/ 4-inch  rotary  feedthroughs.  Viton  O-rings  are  used  for  all 
flanges,  except  for  the  flanges  holding  the  ion  guns,  which  are  of  the  Conflat 
type  with  copper  gaskets. 

2.3  Ion  suns 

The  two  ion  guns  were  obtained  froo  Ion  Tech,  Inc.  (Ft.  Collins, 
Colorado),  Model  2.5-1500-125.  Both  have  2.5ca  diaaeter  grids  and  are  run  by 
2 .5-1500-125EL  power  supplies.  Both  can  be  closed  off  by  shutters  contsining 
ion  probes.  So  far,  only  argon  has  been  used  as  a  sputter  gas  in  either  gun. 

The  gun  aiaed  at  the  target  is  typically  run  at  a  beaa  voltage  of  1500V 
and  a  current  of  70-100nA.  A  feedback  systea  (MKS  250B)  that  regulates  the 
gas  supply  to  the  guns  keeps  the  beaa  current  constant.  The  chaaber  pressure 
during  the  run  is  between  5xl0~5  and  2xl0-*  torr. 

The  gun  aiaed  at  the  substrate  is  used  aainly  for  sputter  cleaning  prior 
to  deposition,  usually  at  a  beaa  voltage  of  300V  and  a  current  of  30mA.  The 
gas  flow  to  this  gun  is  regulated  by  a  three  channel  controller  (MKS  254)  that 
makes  it  possible  to  inject  a  mixture  of  three  gases,  which  could  be  used  for 


controlled  reaction  with  the  spattered  species  at  substrate. 


2.4  Vacuum  system 

The  system  is  pumped  by  a  6-inch  diffusion  puap,  separated  frost  the 
chamber  by  a  liquid  nitrogen  filled  cold  trap.  The  pressure  is  read  by  an 
ionisation  gauge  located  jnst  below  the  vacua  post  of  the  chamber.  Without 
baking  out  the  system,  a  base  pressure  of  3x10  torr  is  obtained. 

2 .5  Deposition  control  system 

Figure  2  shows  a  schematic  of  the  feedback  system  used  to  produce 
compositionally  modulated  films.  The  sputtering  rate  is  monitored  by  an 
Inficon  crystal  sensor  (vertical  model  XFN7504061)  connected  to  an  TTC 
monitor.  The  latter  is  connected  via  an  RS-232  port,  set  at  1200  baud  (bits 
per  second),  to  an  IBM- PC  with  128k  of  memory  and  one  double- sided/double¬ 
density  disk  drive.  The  PC  program  has  as  input  parameters  the  total 
thicknesses  and  densities  of  the  alternating  layers.  When  the  thickness  of  a 
layer  reaches  the  preset  value,  the  PC  sends  a  set  of  pulses  to  a  stepping 
motor  translator  (Superior  Electric  Slo-Syn  type  ST103),  which  makes  the 
stepping  motor  (Slo-Syn  MD93-FC11,  geared  down  by  8.5  using  a  toothed  belt) 
turn  the  target  assembly  by  90°  to  the  other  target.  This  rotation  takes  0.7 
seconds,  and  requires  no  ramping  of  the  speed.  The  instantaneous  deposition 
rates,  total  sample  thi.xness.  deposition  time,  and  the  deposition  rates 
averaged  over  a  preset  time  interval,  are  continuously  monitored  and 
calculated,  and  are  periodically  printed  out  on  a  CRT  monitor  and  a  printer. 

The  PC  has  also  been  programmed  to  alternate  targets  based  on  preset  time 

o 

intervals;  this  is  useful  for  very  short  modulation  repeat  lengths  (1-2A),  for 
which  the  resolution  of  the  crystal  monitor  is  insufficient. 
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3.  EXPERIMENTS 

3 .1  Depoiition  rates 

A  typical  deposition  rate,  as  measured  at  the  sensor  position,  for  pure 

copper,  sputtered  at  a  beam  voltage  and  current  of  1S00V  and  70aA, 

—6  0 

respectively,  and  a  chamber  pressure  of  SzlO  torr,  it  2.5A/sec.  The 
deposition  rate  at  the  substrate  position,  which  could  be  aeasured  quite 
accurately  from  the  z-ray  aodulation  peak  of  the  aodulated  filas,  was  found  to 
be  between  10%  and  30%  greater  than  at  the  sensor  position.  The  thickness 
variation  across  the  substrate  was  also  investigated  by  the  z-ray  technique 
and  found  to  be  less  than  1%/cn  at  the  center  of  the  substrate. 

3 .2  Deposition  of  hoaoteneous  filas 

Binary  aetal  alloys  were  deposited  from  cob posit  targets.  An  arrangement 
consisting  of  horizontal  strips  (1/4  to  1/2-inch  wide)  of  one  aetallic 
component  stretched  across  a  4-inch  square  backing  plate  of  the  other 
component  was  found  to  give  the  aost  reproducible  results.  Electron 
aicroprobe  analysis  showed  the  composition  across  the  2-inch  substrate  to  be 
unifora. 

Amorphous  aetal  films  of  composition  Ni^Zr^,  Cttg2Zrig  and  Cu6JZr35  have 
been  produced.  CnpoZrxo  *nd  ,(r(  aicrocrystalline. 

Amorphous  Si  and  Ge  filas  were  produced  froa  4-inch  square  cast 
poly crystal line  targets  that  were,  respectively,  0.5aa  and  2aa  thick.  Best 
results  were  obtained  by  glueing  these  targets  to  an  aliaina  backing  plate 
with  a  silver-epoxy  conducting  glue  (Ablestik).  The  conductivity  of  the 
targets  was  high  enough  to  aake  neutralization  of  the  bean  was  unnecessary. 

The  argon  incorporation,  measured  by  electron  aicroprobe  analysis,  was  1% 


or  lass  in  aost  filas.  Contamination  by  stainless  steel  sputtered  fron  the 


shield  around  the  target  asseably  could  be  avoided  by  keeping  the  accelerator 
voltage  around  100V,  which  limits  the  bean  divergence. 

The  structure  of  the  filas  was  determined  by  a  combination  of  z-ray 
diffraction,  electron  diffraction,  and  differential  scanning  caloriaetry. 

3 .3  Coanoait ionallv  modulated  filas 

3.3.1  Metallic  filas.  Filas  of  Cu/Al,  Cu/Ni,  and  Cu/Zr  with  repeat 

lengths  between  10  and  lOoS  were  produced  fron  pure  targets  on  glass  or  Si 

substrates.  Tn  aost  of  the  filas  the  thicknesses  of  the  two  types  of 

alternating  layers  were  equal.  The  Cu/Al  and  Cu/Ni  filas,  and  the  Cu/Zr  filas 

o 

with  repeat  lengths  above  6SA  were  polycrystalline.  The  Cu/Zr  filas  with 

o 

repeat  lengths  less  than  65A  were  aaorphous;  which  is  a  aanifestation  of  the 
formation  of  the  aaorphous  phase  by  interdiffusion  in  the  solid  state  between 
the  pure  elements  observed,  for  ezaaple,  in  the  La-Au,  Ni-Zr  and  Ni-Nb  systems 
£11-13]  It  is  also  an  indication  of  a  substantial  aaount  of  interdiffusion 
during  the  deposition  process.  Lin  et  al .  [14]  noticed  a  similar  effect  in 
the  D.C.  sputter  deposition  of  Nb/Co  filas.  Nevertheless,  the  aaorphous  Cu/Zr 
filas  were  still  coa positionally  nodulated,  as  indicsted  by  the  presence  of 
the  characteristic  z-ray  peak.  As  aade,  the  Cu/Al  filas  showed  a  very  strong 
first- ardor  z-ray  (000)  aodulation  peak;  the  weakness  of  the  second-order  peek 
and  absence  of  higher  order  peaks  is  further  evidence  of  substantial  rounding 
»£  the  composition  profile  due  to  interdiffusion  during  deposition  After 
five  months  of  storsge  at  roan  teaperature,  the  intensity  of  the  aodulation 
peak  in  the  Cu/Al  filas  had  decreased  to  a  tenth  of  its  original  value;  this 
I. s  t:>  be  ezpeeted  for  polycrystalline  non-textured  layers,  in  which  nixing  esn 
he  enhanced  substantially  at  low  teaperatures  by  grain  boundary  diffusion 
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along  soring  boundaries  [7] . 

3.3.2  Metal-senicondoctor  films.  Two  sets  of  Wb/Si  films  were  produced. 
Their  average  compositions  and  repeat  lengths  were,  respectively:  Nb^Si^j, 
32X;  and  N^g^Si^,  48A.  The  composition  was  adjusted  by  changing  the  relative 
thicknesses  of  the  two  layers.  This  method  allows  much  better  control  of  the 
average  composition  than  spnttoring  from  a  composite  target.  Our  technique  is 
therefore  quite  useful  for  producing  starting  materials  for  laser  mixing. 

The  Nb  in  the  as-deposited  films  is  polycrystalline  with  very  small 
grains;  the  Si  is  amorphous. 

3.3.2  Semiconductor  films.  Amorphous  Si/ amorphous  Re  films,  with  equal 

layer  thicknesses  for  both  components  were  produced  by  sputtering  from  pure 

o 

polycrystalline  targets,  with  repeat  lengths  of  26,  47  and  64A.  They  all 
shoved  strong  first  order  modulation  peaks.  The  full  width  at  half  maximtm  of 
this  peak  for  the  64X  film,  for  example,  was  0.25°  (using  Cr  Ka  radiation). 
This  corresponds  to  a  10%  variation  in  the  layer  thickness  over  the  course  of 
the  deposition.  The  second  and  third  order  peaks  in  this  sample  were, 
respectively,  20  and  50  times  less  intense  than  the  first  peak;  this  again 
indicates  substantial  rounding  of  the  profile  due  to  interdiffusion  during 
deposition.  Nevertheless,  these  samples  are  very  well  suited  for  the  study  of 
interdiffusion  in  the  a-Si/s-Ge  system  by  further  annealing.  The  sharpness  of 
the  composition  profile  can  be  improved  by  lowering  the  deposition  rate;  the 
width  of  the  modulation  peak  can  be  decreased  by  improving  the  stability  of 
the  crystal  sensor. 

4.  CONCLUSION 

The  system  described  in  this  paper  provides  a  new  method  for  producing 


compositional ly  modulated  films.  Its  advantages  are: 

(i)  simple  and  relatively  inexpensive  design; 

(ii)  flexibility,  in  that  it  allows  deposition  of  metals,  semiconductors 
and  insulators,  without  changes  in  the  system;  it  also  serves  to  deposit  alloy 
films  of  good  uniformity  from  inexpensive  composite  targets; 

(iii)  convenience  of  control,  since  the  crystal  oscillator  feedback 
system  makes  it  possible  to  ''dial  in'*  the  required  layer  thicknesses. 

Section  3  was  intended  as  a  survey  of  the  capabilities  of  the  system. 
Further  details  on  the  characterization  and  the  physical  properties  of  the 
materials  discussed  will  be  presented  in  separate  papers. 
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FIGURE  CAPTIONS 


Figure  1.  Schematic  diagroa  of  the  ion  beta  spattering  cbsaber,  viewed  from 
the  top. 

Figure  2.  Feedback  system  for  the  deposition  of  cob positionally  moduli ted 
materials. 

Figure  3.  Photograph  of  the  ion  beaa  sputtering  systea  (left:  power  supplies 
center:  choaber  and  vacuua  systea;  right:  crystal  monitor, 
computer,  printer,  monitor  translator). 


Figure  2.  Feedback  system  for  the  deposition  of  compositlonally  modulated  materials. 
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